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PREFACE 


T t ih iiiiw over tlurtiy yoiUH siiuw IIi'aviHido'f! opoiational inothoda of 
HolvinK t.lu‘ (liHoit'tit.ial oiinations of pliymc'i woio (irHi publiahod, but 
lut.h(nt<> tiu'y liavo u'cinvocl vory httlo attoutioii fioiu tiiatliotnatical 
phyhU‘iHt« in fjouoial 'I'ho (diiof roaHoii foi tluH huH, I thnik, in tlio lack 
of a foimootod ai-comit, of tho iiiotliodH IIoavimdo’H own work w not 
Ryatomatioally atraiiKod, and m plaooH itn moanniK w not vory oloar 
Bromwioh’a diwuHHion of hw method by moaiw of tho theory of funchiotm 
of a eomph'x vanahio eHtablwhed itn validity , and as a niattor of ])ractu!al 
oonvf'iuoueo tlioio can be little doubt that tho operational motliod is fai 
tho boat for doalinff with the oIuhr of problonm ooncornod It is ofion said 
that it will solve no problom that oannot bo solved ot.lu'rwiso Whothoi 
this is tnio would bo difficult hi say, but it is oertam that in a voiy 
lar^?o class of eases the opoiat.ional motliod will give tho answer in a 
paj^o when ordinal y uiethods take five jiayji's, and also that it gives the 
collect answer when tho oidiiiaiy im'thods, through human fallibility, 
are liable to give a wimig one In partieiilar, when wo discuss tho small 
oHcillatiouH of a dynamieal system with w degrees of freedom by tho 
method of normal eoordinates, wo obtain a detorminantal oipiation of 
tho «th degree to givo tho speeds of tho normal modes To find tho ratios 
of the amplitudes wc must then oomploto the solution for each mode If 
we want tho actual motion duo to a given initial (hsturhanco wo must 
solve a further family of 2 m simultenooiis oi] nations, iinloss spooial 
Himphfyuig cireiimstancos aro present. In tho oporaiional mothod a 
fnimal'oiieiational solution is obtained with tho same ainounbof trouble 
as is neediHl t,o give the petiod equation In tho ordinary mothod, and 
from this the complet.e solution is ohbainahlo at. once by a genoral rule 
of interpndation For c(mtinu(mH systems tho advantage of tho opora- 
tiimal mi'thod is evmi greater, for it gives hetli iieriods and amplitudos 
easily in prohleuis where the amplitudes cannot ho found by the ordmaiy 
mothod without a kuowlodgo of some tlieorcun of expansion in normal 
functions analogous to b'enrier’s theorem. Hoat conduction is also 
especially convomently treated by fiporataonal methods. 

Smeo Bromwich’s discussion it lias often been said that tho operational 
mothod is only a shorthand way of writing contour integrals It may bo ; 
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but at least one may reply that a shorthand that avoids the noceshity 
of writing i I dKin every line of the work is worth while Oonneciocl 

27nJ 0— too 

With the saving of writing, and perhaps largely because of it, is the 
that the operational mode of attack seems much the more natural when 
one has any familiarity with it After all, the use of contour lutegralh 
in this connexion was introduced by Bromwich, who has repeatedly dc^ 
dared that the direct operational method of solution is tho better of 
the two 

My own reason for writing the present work is mainly that J liiwo 
found Heaviside’s methods useful in papers already publislicd, and shall 
probably do so again soon, and think that an accessible account of them 
may be equally useful to others In one respect I must offer an apology 
to the reader Heaviside developed his methods mostly in relation tt> tins 
theory of electromagnetic waves Having myself no qualification to writhe 
about electromagnetic waves I have refrained from doing so , but m t.ho 
operators occurring in the theory of these waves ate mostly of typOH 
treated here I think the loss will not be serious It can m any (‘.ano ho 
remedied by reading Heaviside’s works or some of the papers'll! tho h«i 
at the end of this tract 

A chapter on dispersion has been included The operational solution 
can be translated instantly into a complex integral adapted for evalua- 
tion by the method of steepest descents, a short account of tho latter 
method has also been given, because it is not at piesent very accoHHible, 
and IS often incorrectly believed to be more difficult than the method of 
stationary phase Two cases wheie the Kelvin first approximation to the 
wave form breaks down are also discussed 
My indebtedness to the writings of Dr Bromwich is evident from tlu^ 
references in the text In addition, the problems of 4 4 and 4 5 are taken 
directly from his lecture notes, and several others are included largely 
as a result of conversations with him 
My thanks are also due to the staff of the Cambridge UmvorHiiy 
Press for their care and consideration during publication 

HAEOLD JEPKEIYB 


St John’s Collbob, 
Cambridub 
1927 Jul^ 19 



(M)NTENTH 


Prftiic*t^ . , . * page v 

, . * . , vh 

Clifip I* Nntiuim . . 1 

li (lomplttx TluKiry ... . * 19 

liL I9iy«iml ApplicatiouH; One ludepeiidout Variable 27 

IV, Wava Motion in One DiuienHuni , , 40 

V. (Jeiidiinfion <d’ Heat iu One Dimonaiou 54 

VI frobleiim with Hpbeaeal (^r Oylmdneal Symmetry , 06 

VII l)iK|)im4mi , . 75 

VIII, Be^^e! BtuirtieiiH . , , . , 86 

Note On the Ncttalimi fer the Emm Fuimiion or Pro- 
bability Iiitei^ral . 94 

fiiteijiretatioim c»f the Pniieipal Operatora . . 96 

BiMifigraphy . , , . , , 97 

Index k) Autlnmn . , . . , 100 

Suhjirt Index 


101 




CHAPTER I 

FUNDAMENTAL NOTIONS 


1 . 1 . Let us consider the linear differential equation of the first order 




( 1 ) 


^here B and 8 are known functions of bounded and mtegrable when 
0 ^ ^ Suppose further that y = yQ when x-Q Let Q denote the 
operation of integrating with regard to x from 0 to x^ so that 




( 2 ) 


Perform the operation Q on both sides of the eq[uation (l) Then we 
find 


( 3 ) 

and the right side vanishes with x This can be rewritten in the forms 
{l-QB)y^y(,‘¥Q 8 , ( 4 ) 

y=y,+ QS+QBy ( 5 ) 

These are both equivalent to the original differential equation (1) to- 
gether with the given terminal condition When we write QRy we mean, 
of course, that B is to be multiplied into y and the product integrated 
with regard to x from 0 to ^ But the whole expression for y may be 
substituted in the last term of ( 5 ), giving in succession 

y-^yo'^QS-^QB (3/0 + QRy) 

=3^0 + QS + QB (3/0 + QS) + QBQB (3/0 + Q 8 + QRy) 

= (yo + QS) + QR (yo + QSO + (yo + QS) 

+ QRQBQR (3/0 + Q 8 ) + (6) 

on repeating the substitution indefinitely We have to show that the 
infinite series (6) converges and that it is the correct solution In 
evaluating each term R is supposed to be multiplied into the whole 
expression after it, and Q to operate on the whole expression after it 
Suppose that within the range considered 

\B\^A, \yo+Q 8 \^B, . ( 7 ) 


where A and B are finite Then the absolute value of the second term 


j 


I 
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IS less than ABoc, that of the third than that of the fourth 

than A^Bix^lB ^ , while the general term is less than A^Bafln ’ The series 
therefore converges at least as fast as the power series for exp Aa^ It 
therefore represents a definite function, and on substituting it in (5) we 
see that the equation is satisfied for all values of a) ^ also the solution 
reduces to % when w = 0y as it should Thus (6) is the correct solution. 
The solution can also be written 

y^(l-^QB + QBQB^QBQBQB+ ) (yo + QS) ( 8 ) 

The operator between the first pair of brackets is the binomial expansion 
of (1 - carried out as if QB was merely a number Since yo + 

is a determinate function, we can write the solution in the form 

3 / =■ iZTqB 

provided that yo + QS is evaluated first and that the operator (1 - QBY^ 
is expanded by the binomial theorem before interpretation In fact (9) 
IS merely a shorthand rule for writing (8) But on returning to (4) we 
see that (9) is also the solution of (4) earned out as if 1 - QB was a mere 
number 

It IS evident from (8) that the values of 8 for negative values of x do 
not affect the solution provided yo is kept the same Suppose then that 
8 was zero for aU negative values of and that y was zero when x=Q 
The solution would be 

y^(l + QB+QBQR+ )Q8, . (10) 

and this solution would be unaltered if the lower limits of all the 
integrals were replaced by - oo But if now we add to ^ a constant yo/f 
for aU values of x between 0 and i, Q8 will be increased by yo for all 
values of x greater than and (10) will be converted into (8) If then $ 
tends to zero, yo remaining the same, y^ will tend to yo, and we recover 
the solution (8) with the original initial condition The physical interest 
of this result is that it corresponds to our notions of causality Suppose 
that the independent variable x is the time, and that y represents the 
departure of some variable from its equilibrium value , then B represents 
a property of the system and 8 an external disturbing influence If the 
system was originally in its equilibrium state, the form (10) exhibits 
the disturbance produced by the external influence after it enters If it 
was undisturbed up to tune zero, the part of (8) dependmg on 8 repre- 
sents the effect of the fiinite disturbances acting at subsequent times, 
while the part depending on yo represents the effect of the impulsive 
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disturbance at time 0 required to change y suddenly from zero to If 
we like we can separate the solution into two parts 

{i-qRY^y, and (i-qBT^qs, 

and say that the first represents the effect of the initial conditions and 
the second that of the subsequent disturbances 


1-11. The method just given can be extended easily to cover many 
equations of the second and higher orders Thus if our equation is 


dx^ 




( 1 ) 


with y-ya and ^ when x = 0, we find by integration 


= + ( 2 ) 

y-y,>-xy^ = QQ{By + S) ( 3 ) 

This leads to the solution in a series 

y = {l + q^It + + ) (y„ + xy, + (4) 

'wliere Q^f^x) = jV(#) didt (5) 


1.12. As a special example, consider the equation 

with y = l when x = Q Then carrjang out the process of 1 1 (6) we get 
3r = (l + aQ + a=*Q=*+ )l (2) 

( 3 ) 

the ordinary expansion of exp ax Or take the equation 
with y -I and dyjdx = 0 when x = 0 We can infer 


1 + ax + + 


2 /=!- ( 5 ) 



the ordinary expansion of Jq (x) 
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1 . 2 . The foregoing method is due originally to J Caq[u4*^ it is a 
valuable practical method of obtaining numerical solutions of linear 
differential equations Its extension to equations of any order, or to 
families of simultaneous equations of the first order, is more difficult 
unless special simplifications enter But if the equations have constant 
coefficients, which is an extremely common case in physical applications, 
a considerable development is possible This arises from the fact that 
the operator Q obeys the fundamental laws of algebra Thus if a is a 


constant and u Bind v known functions of co, 

= ( 1 ) 

Q(u + 'd) = Qu+ Qv, (2) 

Q^Q^u = Q^Q^u = ( 3 ) 

Consequently Q behaves in algebraic transformations just like a number 
If for instance we have two operators /(Q) and g(Q) both expressible 
as sums of integral powers of Q, thus 

^(Q) = ao + aiQ + a2Q^ + a^^+ (4) 

9{Q) = bo’^biQ + b2Q^ +bsQ^+ ( 5 ) 


where the a’s and 6's are constants, let us for a moment replace Q by a 
number z small enough to make the series converge absolutely Form 
the product series 

/(^) ^ (^) = (^0 + + ^ 2 ^;^ + ) {bo + biZ‘¥b2Z^+ ) 

== Co + + ^2^^ + (6) 

say Consider first the case where the series are both polynomials Then 
if 8 IS an integrable function of ^ 

/(Q) g(Q)S^(co^G,Q-^c2Q^+ )S (7) 

For the left side means 

{ao + + asQ" + ) (bo8+b, Q8+ b^Q^S^ ) 

= ^0 (boS+b^QS {bo Q8 + b^ Q^8 + ) 

+ {boQ^S + ) + (8) 

* Ltouvtlle^s Journal, (2), 9 (1864), 185-222 Further developments are given by 
Fuchs, Ann d Matem , (2), 4 (1870), 36-49 , Peano, Math Ann , 32 (1888), 450-456 , 
H F Baker, Proc Lond Math Soc , (1), 34 (1902), 347--360, (1), 35 (1902), 333- 
378, (2), 2 (1905), 293-296, Phil Tram , A, 216 (1916), 129-186 Caqu^ considera 
only a single differential equation, hut notices that the operator in the result is 
in the form of a binomial expansion For other operational methods based on these 
principles, but applicable to equations of order higher than the first, or to families 
of equations of the first order, the above papers of Prof Baker may be consulted 
Physical apphcations are given by W L Cowley and H Levy, Phil Mag , 41 
(1921), 584-607, Jeffreys, Proc Lond Math Soc , (2), 23 (1924), 454 and 465,. 
M NBA S , Geoph Suppl 1 (1926), 380-383 
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using (1) and (2), and then using (3) we can collect the terms involving 
the same power of Q and obtain 

ci^oboS + (a^ hi + ai QS + (<350^2 + (iihi + 5 o) + ( 9 ) 

which IS by definition the same as 

(co + Ci^ + C2^"+ )8 (10) 

When the series are infinite, their convergence may be established 
easily If the senes for/(; 2 ;) converges absolutely when \z\^r, then a 
number M must exist such that ^ M for all positive integral values 
of n Thus 

( 11 ) 

and also if for all values of | /S' | ^ (7, where (7 is a constant, 

«•««»!; (12) 


Hence It /{«)«- i «.«•«, (13) 

71=0 ^ 

each term is less than the corresponding term of the series 



(U) 


which converges for all values of x however large So long as f{z) is 
expansible about the origin in a convergent power series, however 
small its radius of convergence may be, the expression /{Q) 8 will be 
the sum of an absolutely convergent series however great x may be If 
f(z) and g{z) are both expansible within some circle, their product 
series will also converge within this circle and the expression/(§)^ (Q) 8 
will be an absolutely convergent series however great x may be 
Provided with this result we can now easily extend the result (7) to the 
case where /( 2 ;) and g{z) are infinite series, by methods analogous to 
those used to justify the multiplication of two absolutely convergent 
power series 


1.3. We can now extend these methods to the solution of a family of 
n simultaneous differential equations of the first order with constant 
coefficients Suppose the equations are 

Vi + ^12 ^2 + + yn = 8i 

^2iyi + 622 y 2 + +^2wy7l=/S'2 

y2 + + BnnVn = 8^ 
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where the y’s are dependent yariables, x the independent variable, era 

denotes ars-r- + Ks, where ars and hrs are constants, and the are 
ax 

known functions oi x We do not assume that 

Chrs ~ ei<rr , ~ j 

but we do assume that the determinant formed by the a's is not zero 
When x-0^ 2 /i = 'Wi, and so on, where the are known constants 
First perform the operation Q on both sides of each equation We have 

Q,er&y8 ~ f f^rs "f" hray^^^ 


~ fvB ys ^ra 3 ( 2 ) 

where fra denotes the operator + Then the equations and the 
initial conditions are together equivalent to the equations 
fii yi +^12 y^ + ’^fin yn — '^i'^ 


f niyx'^fn'iy%'^ frmyn — 

where 'yr = ^rl«^l+«r 2 «^ 2 + -^amUn (4) 

The general equation can be written compactly 

^8 fra ~ '^'T Q,^t ( 5 ) 

Now let D denote the operational determinant formed by the /’s, 
namely, 

fil fi2 fin 

fi ^22 fn , . 

(6) 

fnl fn% fnn 

If this determinant is expanded by the ordinary rules of algebra and 
equal powers of Q collected, we shall obtain a polynomial in Q The 
term independent of Q is simply the determinant formed by the a’s, 
which by hypothesis does not vanish Now let J^rs denote the minor of 
fra in this determinant, taken with its proper sign JF'ra also is a poly- 
nomial in Q 

Now operate on the first of (3) with Fjs, on the second with jPL, and 
so on, and add Then in the sum the operator acting on y^, say, is 

^rFrafrm (l) 

lim = s, this sum is the determinant i) , if m = 1 = 5 , it is a determinant with 
two columns equal and therefore is zero The resulting equation is 
therefore 


I)ya==^rFra(Vr’^QSr) 


( 8 ) 
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Now if all tte 8’^ are bounded and mtegrable within the range of values 
of o) contemplated, the expression on the right of (8) is also a bounded 
mtegrable function of Also since the function 1) (z), obtained by 
replacing ^ in D by a number z, is regular and not zero at z = 0, the 
function II D (z) is expressible as a power series in ^ with a finite radius 
of convergence Define as the power series in Q obtained by putting Q 
for in the series for 1/D (z) Then operate on both sides of (8) with D~^ 
We have 

= B-^XFrs (Vr + QSr) ( 9 ) 

But, since series of positive integral powers of Q can be multiplied 
according to the rules of algebra, D“^D gives simply unity, and we 
have the solution 

ys = D'^'^XFrs + QSr) ( 10 ) 


This gives a complete formal solution of the problem Its form is often 
convenient for actual computation, especially for small values of x, but 
it can also be expressed in finite terms D is, as we have seen, a poly- 
nomial in Q of degree n at most, while Frs is a polynomial m Q of 
degree - 1 at most Our solution is therefore of the form 


y OS 

D(Q)^ 


( 11 ) 


where and if/ are pol 5 niomials whose degree is ordinarily one less than 
that of D Since the determinant formed by the a’s is not zero we may 
denote it by A, and then D will be the product of n linear factors, 
thus 

D = A (l-aiQ)(l-a2Q) (l-a^Q), (12) 


where the as will ordinarily be all different Then </>(Q)/D(Q) can be 
expressed as the sum of a number of partial fractions of the form 


L 

1-aQ 


( 13 ) 


But by 1 12 (3) this is the same as The part of the solution arising 
from the ^'s can therefore be expressed as a linear combination of 
exponentials 

The justification of the decomposition of <t>(Q)/D(Q) into partial 
fractions is that this decomposition is a purely algebraic process Hence 
if the partial fractions and the original operator are all expanded in 
positive powers of Q, and like powers of Q are collected, the coefficients 
of a given power of Q will be the same in both expressions, and the 
eq.uivalence is complete 



® fundamental notions 

Exceptional cases will occur if D contains no term in or if two 
or more of the as are equal If the term in is absent the expansion 
of the operator inj)artial fractions will usually contain a constant term 
ft® term in § ^ is also absent we must divide out, and the expan- 
sion in partial fractions will contain a term in Q, this will give a term 
in X on interpretation 

If several of the a’s are equal, the expression in partial fractions 
mil involve terms of the form M(l-o,q)-- These can be interpreted 

y irect expansion , but another method is more convenient Starting 
with 


let us differentiate r—l times with regard to a 


(l-a^r 





We find 


(14) 

(15) 


so that 


_^-i _ ^-1 

(l-aQJ (r-l)! 




( 16 ) 


A rather different form of resolution into partial fractions from the 
ordinary one is therefore necessary if each fraction is to give a single 
term in the solution Instead of having constants m all the numerators 
we must have powers of Q, the power of Q needed in any fraction being 
one less than the degree of the denominator But if we write for Q 
the fraction in (16) is algebraically equivalent to pjijp - a)*- In tbs form 
the power of p in the numerator is independent of the degree of the 
denominator, and it is easier to resolve into partial fractions of tbs form 
tLan mto tliosG iiiYolviiig Q directly 

The above remarks apply to the interpretation of the effect of the 
initial conditions, that is, of the first term on the right of (11) To find 
the effect of the terms in the ^’s we may expand the operators simi- 
larly To interpret (1 - aQ)-^ Q8 m finite terms, we note that it is the 

solution oi -£^-ay=8 that vanishes with x Tbs is easily found by the 


ordinary method to be 


y = ^Q(8e-'^), (17) 

which IS the bterpretation required 
To sum up, we can solve a family of equations of the type (1) by 
first int^ating each equation once from 0 to a; with regard to x, allow- 
conditions Tbs gives a set of equations of the type 
(3) 1 he subsequent process for deducing the operational solution (10) 
IS exactly the same as if the operators /„ were numbers, and the 
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ordinary rules of algebra were applied The solution can be evaluated 
by expanding the operators in ascending powers of Q and evaluating 
term by term, this in general gives an infinite series Alternatively it 
can be obtained by resolving into partial fractions and interpreting each 
fraction separately, this gives an explicit solution in finite terms 

1.4. Heaviside’s method is equivalent to that just given , it differs 
in using another notation, which is not quite so convenient in a formal 
proof of the theorem, but is rather more convenient for actual applica- 
tion In Heaviside’s notation the operator above called Q is denoted by 
p At present we need not specify the meaning of positive powers of 
P ! negative integral powers are defined by induction, so thatyi“’® denotes 
W 6 have noticed that the passage from 1 3 (3) to 1 3 (10) is a 
purely algebraic process Consequently if all the equations (3) were 
multiplied by constants before the algebraic solution the same answer 
would be obtained Suppose then that we write the general equation 
1 3 (5) in the form 

'^a{ars + b„p~^)y,= \arsUs+p-^Sr . ( 1 ) 

Multiply throughout by y? as if this were a constant We get 

Ss (flraP + ^rs) Pa — '^a^raP'^a "t' (2) 

On solving the n equations of this form we shall obtain a solution 
identical with 1 3 (10) except that p-^ will appear for Q, and both 
numeiator and denominator will be multiplied by the same power of p 
If the operators in the solution are expanded in negative powers of p, 
&nAp-^ IS then interpreted as Q, the result wiU be identical with that 
already given Comparing (2) with the original equations 1 3 (1) we see 
that the new form of our rule is as follows 
Write p for djdx on the left of each equation, to the right of each 
equation add the result of dropping the &’s on the left and leplacing 
the ys by their initial values, solve the resulting equations (2) by 
algebra as if p was a number , and evaluate the result by expanding 
in negative powers of p and interpreting as the operation of 
integrating from 0 to ^ 

This IS Heaviside’s rule In what follows the equations (2) will usually 
be called the subsidiary equations 

To obtain the solution explicitly, we put era for arap + denote our 
determinant 

^12 €in 

^21 ^22 ^2ri 

&n2 0nn 


I 
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l>y A, and denote the minor of e„ in this detemnnant, taken with its 
proper sign, by E^b Then the solution is 


E 


(3) 


Since the determinant formed by the a’s is not zero, A is of degree n in 
p, while Era IS at most of degree w — 1 The operators can therefore be 
expanded in negative powers oip, as we should expect, positive powers 
do not occur All terms after the first vanish with a>, the first is 


S -^T8 


(4) 


where is the minor of ars m A But 


= 0 (m4=5)J ' ^ 

and (4) reduces to as we should expect This verifies that the solu- 
tion satisfies the initial conditions 


1.5. To interpret the operational solution 1 4 (3) in finite terms we 
req^uire rules for interpreting rational functions of^ operating on unity 
and on other functions We have already had the rules 

P~'^=Q, = O' , and so on, (1) 

= 07, p-^1 = Qo 7 = ~ , and in general = 2! ( 2 ) 

If unity IS replaced by Heaviside’s ‘unit function/ here denoted by 
H (jx^^ which IS zero for all negative values of x and 1 for all positive 
values, we shall still have 

= (3) 

when X is positive, hut it will vanish when x is negative We can also 
replace the lower limit of the integrations by — 00 without altering this 
interpretation 

Again, we shall have when x is positive 


( 4 ) 

p _ 

{p - a)" (1 - a§)» ~ (re _ 1 ) I (5) 

^=^-1 = ^- 1 , ( 6 ) 

where the function operated on may be either unity or E{x) In the 
latter case all the operators wiU give zero for negative values (Ax 
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The operators m 1 4 (3) are of the iQicmf^p^lFQp)^ where /{p') and 
F {p) are polynomials in p, and f {p) is of the same or lower degree 
than F{p) If F{p) is of degree n it can be resolved into n linear 
factors of the form p — a Then provided that the a's are all different 
and none of them zero we have the algebraic identity 
f{p) _ /(O) . ^ /(a) 1 


whence = tJ'Ps + 

F(p) F(0) ^^ccF'(a)p-a 


aF' (a) p - 
/(“) p 


( 7 ) 

( 8 ) 


If this operates on unity or H (x) we have therefore for positive values 
of a? 

J^(p) Flo) ^‘^aF' (a) ^ 

To justify this we notice as before that (8) is a purely algebraic iden- 
tity, and therefore if both sides are expanded in negative powers of p, 
beginning with constant terms, the expansions of the two sides will be 
identical, and on interpretation in terms of integrations will give the 
same result 

The formula (9) is usually known as Heaviside’s expansion theorem , 
but as Heaviside’s methods involve two other expansion theorems* it 
will be called the ^ partial-fraction rule ’ in the present work 

If some of the a’s are equal or zero, the expression (7) considered as 
a function of p will have a multiple pole, and its expression in qiartial 
fractions will contain terms of the form (p — a)~®, where s is an integer 
greater than unity, and a may be zero Then f (p)IF {p) will contain 
terms of the forms or pUp — ay, which can be interpreted by 

means of (2) or (5) 

By means of these rules we can evaluate all the expressions for the 
part of ys in 1 4 (3) that depends on the initial values of the y’s If the 
8 s are constants for positive values of x, as they often are, the same 
rules will apply to the part of the solution depending on them If they 
are exponential functions such as the easiest plan is usually to 
rewrite this as p/lp — fP) and reinterpret Thus 


f(pl 

F(p) 




f{p) 

F{p)p-iA. 


( 10 ) 

( 11 ) 


Namely, expansion in powers of Q or jp-i and interpretation term by term, 
and expansion in powers of e-P\ where h is a constant, as m 4 2 
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If 8 IS expressed as a linear combination of exponeniialB wo can apply 
this rule to each separately This is applicable to prac ii<*ally all func- 
tions known to physics 

Alternatively we can resolve the operator acting on 81 into piirtial 
fractions and interpret p'^^S by integration and other fractioiw by tlio 
rule of 1 3 (17) 

-i— S = <r (12) 

p-o. Jo 

This completes our rules for solving a set of linear ec) nations t»f the 
first order with constant coefficients In comparison witli the ortliniiry 
method, we notice that the rules aio direct and lead iniinediatidy to a 
solution involving operators, which can then bo evaluated cnnnpleiely 
by known rules If it happens that we only rcipiuo the variation of one 
unknown explicitly, we need not inter|)rot the Holutions foi the others. 
In the ordinary method we have to find a complementary function and 
a particular integral separately To find the former we aHHUiiie a solu- 
tion of the form and on substituting in the diffenuitml cHiuti- 

tions with the omitted we find an ecpiaiion of (amsiHfem*y to 
determine the n possible values of a, and tlio ratios of the c«»rrco 
spending to each The particular integral is then foiiml by some nud limit 
but it does not as a rule vanish with x The actual valines of the X’s 
are still undetermined, and the value associated with each « must be 
found by substituting m the initial conditions and again solving a set 
of n simultaneous equations The labour of finding the oqnat 4 on for the 
a’s and the ratios of the Vs corresponding to one of them is about, the 
same as that of finding the operational solution lu I IcaviHide’s mcdluHi ; 
the rest of the work is avoided by the opoiational mothod. Further, tf 
some of the a's are equal or zero considerable compheatioiiH are intro 
diiced into the ordinary method, but not into the operational one. 

1 . 6 . The above work is applicable to all cases whore A, the deter 
mmant formed by the a% is not zero Wo can show that if A m zero 
there is some defect m the specification of the system. A system is 
adequately specified if when wo know the values of tlio diqiendcuit 
variables and the external disturbances the rates of change of the 
dependent variables are all determinate, and if the initial values of tlie 
vanables are independent of one another Now if A m zero, UA us 
multiply the r’th equation by Ars, the minor of in A, and mU up 
for all values of r The coefficient of dyS^ m the sum is A, whicdi is 
zero , that of any other derivative is a detoiminant with two i?oliiwns 
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i;{ 


equal, and w alno zero. TIuih tho derivatives ilisappear ontiroly from 
tlio sum, and wo aro left with a relatiion Imtweeii tim t/'h and the *V’h. 
If the coeflioiouts of the »/« also vanish, fhewe must ho a pormaiiont 
relation botweeu the /S"h, and oiio of Hm equations wi* started witii is a 
mere logical consequeiico of the otlu^rs, then wo have not enough 
equations to determine the dorivative.s If the eoeilieients of the i/'a do 
not vanish, wo can put ir 0 and obtain a relation ladwoen the initial 
values, which aro therefore not independent.. 


1.7. The motliod m most easily evh'iided to eqimtions of higher 
ordei by breaking them up into eipiat.ious of the tirst order. 'I’hu.s if we 
have an equation of the second older such as 


f/“// 


+ « 


r/y 


e Ifi/ 


as; 


we introduce a now vanablo z given by 


(fjj 


( 1 , .. 


■(0 


( 2 ) 


and the original ociuation can ho roplacisl by 

•• Cl) 

Wo have now two equations of the first, ordi'c in // and i. If initially 
y —^0 titul di/jdj the siilmidiaiy eiiuatioiiH .III) 


RV « 


H) 

(p s «) ; 1 hi/ pi/, t 


(h) 

Solving by algebra wo find 

(■/!« 1 (l/l) l/„ 1 1 

.S 

• C'O 

,,, , ' o t (ip 1 h 

To interpret, put 



p'^-ilt/nh (p a)(p 

/4 . 

.. f7) 

and apply the pai tial-fraction rule. We find 

•'''"a A'/u)''" «.'Ad 





...(«) 


whioh IS easily shown 
the solution required. 


to sati.tfy all the eoiiditions and thereforo to Isi 



14 


FUNDAMENTAL NOTIONS 


1,71. A few illustrative examples may be given (p is written for 
djdcc from the start ) 

1. 0^ + 4^+ 3) 3^ = 1, yo= 3, 3^1 = - 2 
We consider the subsidiary eq_uation 

+ 4^ + 3) 2/ = 1 ~ 2jt? + 3 + 4p) 

= 3^^ -f lOp + 1, 

3^" + 102? + 1 

y o+i)(j?+3) 

on interpreting by the partial-fraction rule 

2. (p= + 5y) + 6)2/ = 12, yo=2, 3^1 = 0 

Consider (jp'® + 5^ + 6) y = 12 + 2 + bp) 

Cancelling the common factor, 

2/=2 

3. (i> + 3)2/ = 6-“", yo = 0 

This can be written (i> + 3) 2 / = , 

jP + 2 

2' = (p + 2fo + 3) 


+ ?^^) 2/ = 0 

p'^y^i+pyi sinwa? 

y = ■:=y,C 0 S^-fy,— 


5. {p^^yy-=o^e-^, yo=0, 3 ^ 1=0 


The expression on the right is equivalent to the operator 2pl(p + 2)^ 
Hence 


_ 2jt? _ 2i??* 

^“0 + 2)'“'4T 



^q-2sd 


We notice the advantage of Heaviside’s method in avoiding the use 
of simultaneous equations to determine the so-called ‘arbitrary’ constants 
of the ordmary method In particular in the second example the data 
have a property leading to a simple solution Heaviside’s method seizes 


upon this immediately and gives the solution in one line, with the 
ordmary method the simultaneous equations for the constants would 
have to be solved as usual 



PUITDAMENTAL NOTIONS 15 

In general we may say that the more specific the problem the greater 
will be the convenience of the operational method 

1.72. The method just given can be extended easily to a family of 
equations of the second order If the typical equation is again written 

= ( 1 ) 
d? , d 

( 2 ) 




where now 
we can introduce n new variables Zi, z^, 


Zn given by 


da> 


Zr = 0, 


( 3 ) 


thus treating the first derivatives of the y’s as a set of new variables 
Ihen the equations (1) are equivalent to 

dz 


and (3) and (4) constitute a set of 2n equations of the first order Suppose 
also that when a! = 0, ys = u„ Zs = v, According to our rule of 1 5 we 
must replace djdx by p and add pur to the right of (3), and 
to the right of (4) We have now to solve by algebra, and may begin by 
writing the revised form of (3) “ ^ 

^r=p(yr-Ur) ( 5 ) 

When we substitute in the modified form of (4) we get 

^^<^’-sP^Q/3--Ua)+'S„brsp(Sa-U,) + -^,Crs^a = Sr + XaraPV„ ( 6 ) 
or, on rearranging, 

2, (araP^ + braP + Cra) Va = %a {a„f + hrap) M. + 2, a, a pVa + Sr (7) 

W-e have thus ra equations for the y’s to solve by algebra, and the 
solution can be interpreted by the rules 

In Bromwich’s paper ‘Normal Cooidinates in Dynamical Systems’ a 
method equivalent to the operational one is applied directly to a set of 
second order equations First order equations, however, arise in some 
problems of physical interest, and merit a direct discussion This as we 
have seen, is easily geneiahzed to equations of the second order ' 

1.73. In the discussion of the oscillations of stable dynamical systems 
by the ordinary method of normal coordinates there is a difficulty when 
the determinantal equation for the periods has equal roots This does 
present method If the system is not dissipative, 
and the determinant formed by the .’s as defined in (2) has a multiple 
factor {p + «-)’•. we know from the theory of determinants that every 
rst minor has a factor (p^ + a®)’ -i On evaluating the operational solution 
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of (7) we can therefore cancel the factor (j?^ + from the numerator 
and the denominator, and we are left with a single factor {p^ + a^) in the 
denominator Thus the part of the solution depending on the initial 
conditions is of the trigonometric foim as usual terms of the forms 
t cos at and t sin at do not arise 

1.8. In all the problems considered so far the operators that occur in 
the solutions are expansible in positive powers of Q, or in negative powers 
of and we have seen that so long as the series obtained by replacing 
Q by a number have a finite radius of convergence, however small, the 
operational solution is intelligible in terms of the definitions we have 
had But we have not defined p as such, because we have only needed 
to define its negative powers, and p cannot be expanded m terms of its 
own negative powers Then has p any meaning of its own^ Since it 
replaces djdx when we form the subsidiary equation we may naturally 
suppose that p means djdx, and this is the meaning sometimes actually 
attributed to it, but care is needed We recall that when the subsidiary 
equation is formed a term like ^ 3^0 appears on the right, but dyojdx is 
zero, so that if we pushed this interpretation too far we should be faced 
with the alarming result that the solutions of the equations do not 
depend on their initial values The fact is that though the operators 
djdx and Q both satisfy the laws of algebra and are freely commutative 
with constants, they are not as a rule commutative with each other 
Thus /7 /7 /*« 

Q-^f{x) = j^f\a;)dx=f(x)-f{0) ( 2 ) 

Thus the operators p and Q are commutative if, and only if, the function 
operated on vanishes with x It appears from (1) that djdx undoes the 
operation Q if djdx acts after Q With this convention we can identify 
p, the inverse of Q, with djdx When p and Q both occur in an operator, 
the Q operations must be carried out before the differentiations^ 

This result explains why some of our interpretations differ from those 
given in text-books of differential equations for the deteimination of 
particular integrals For instance, we have interpreted — a) as 

” - 1) In the ordinary method, due to Boole, we expand this operator 

in ascending powers of^, and get 



* Cf Heaviside, ElectroTnagnetic Theory, 2, 298 
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This assumes that the sum of all the terms of the senes after the first 
IS zero Jout 

Vi? -a aj^ a(p^a)^ W 

If the operator on the right is wntten i ^ ^ 1 , it is clearly zero. 

But if It IS written ^ j? 1 , it gives by our rules ^ The differ- 
ence anses from the fact that the operators p and (p - a)-^ are non- 
commutative y “io nun 

We may notice, incidentally, that whereas the senes in powers of 0 
give convergent series on interpretation, the same is not true of the 
corresponding series in powers of p For instance, if /S' = (a? - 1 )’-, where r 

IS firactional, the senes for S m ascending powers of p diverges 

hke Apart from the greater internal consistency of 

hive apphcation 

=/(«;) + A/' (a!) + ^/"(^) + 

—f + A), 

S'® IS useful in enabling us to find a theorem that plays 

n Heaviside s method a part closely analogous to that played by Fourier’s 
integral theorem in the ordinary method Instead of tngonometric 
fonctions, Heaviside treats as fundamental the function here called 
HW, which IS zero for negative, and unity for positive, values of a;* 
Ihis function being discontinuous, we cannot at once apply Taylor’s 

we 1 ket and the results obtained by applying (5) may be regarded as 
the limits of results proved for these Then we shall write 

= 1 when oo-h^ 

= 0 when 




( 6 ) 


J 


* Proc Roy Soc , A, 62, 513 

t E g J (ErfXs + l) where \ may be made indefinitely great 
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A fuller justification of this step will be given in the next chapter Then 
if 

H (x)=l if < ^ 2 1 

~0\ix<hi I (7) 

and = Q li x>7i2 J 

If further \ <li^^ 

[{er'^K 4 - 

= 0 for for hT,<cc<h^, and so on (8) 

If then the subdivisions of the interval to become indefinitely 
numerous, we can approximate as closely as we like to a function /(^) 
by a sum 2/ (^) - e~^^) H (x) N ow with a further proviso about 

replacing H{x) by a continuous function and later proceeding to a 
limit, we can replace this sum by an integral Then 

/(^) = r f(J^)d = r (9) 

J h~ 00 7-00 

If the integration with regard to h is carried out, we obtain a function 
of p operating on H {w) In this way a given function of x can be 
expressed in the operational form^ 

We notice that if Ji is positive 


={: 


and 


p lx when x>0 

0 when x<k 
x-h when x>k, 

i e~p^H{x) H(x-K)= j H{x-}i)dx 

= 0 when x<h, 

— x—h when x^h 

Thus the operators ~ and can be commuted provided h is positive 

P 

IrTf /Owhen^<-A 

p liT + Awhen«;>-7^, 

0 when x<0 
.X when ^ > 0 

Fortunately symbolic solutions seldom contain operators of the form 
where Ji is positive, so that the fact that this operator is not com- 
mutative with - gives little trouble in practice 
P 


But 


and 




Heaviside, Electromagnetic Theory, 3, 327 
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CHAPTEE II 

COMPLEX THEORY 

‘"‘ft""””® '»»“i ®» fe the reenltt of opomtaa 

rferi fU^S °T '‘“‘**5' 

rules, as lollows If </> is an analytic function, 

( 1 ) 

K ' ' (2) 

n t e ormer definition the integration extends around a large circle 
m to cmpta ptao lo to Uto to mtogeotoo „ aloogl o^o 
om c 1 00 to c + i CO, where c is positive and finite, such that all 
singularities of the integrand are on the left side of the path (that is 
the side including - oo ) These rules are both due to Bromwich* 
Considering first (1), we have 

^ ~2?t/o^i^'^> (3) 

integrand is at the origin, where the residue 
s ^ Hence when ^ is a positive integer 






£- 


(4) 


Also we see that 

^ce to to »„ m to fimte part of to p£., 

'tte mt«pr,t.t.on I 5(5) otp/^ll). gko follows m.mod.oWy frL 

W. "a- 

/(k) 

^(p) ‘2^tJc>c 1 (AC (K-a„)' 


^dK 


= ^ + 2 /(«) 


(7) 

( 8 ) 


on evaluating the residues at the poles 0, a^, a^, 

4 o!- 28 '“'' Systems,- Pmo Lafl "Matk Sac , 15 (1916). 
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We can also prove that 1 4 (3), ■when interpreted by the rule (1), is 
the correct solution of eq,uations 1 3 (1) If we denote by Crs («), A (k) 
and JSr, (ic) the results of replacing p by « in e™, A, and the part 
of 1 4 (3) not involving the &s is eq[uivalent to 


Substituting in the dififerential equations 1 3 (1) we hnd that the left 
side of the m’th equation is 




(10) 

But 

5s Prs (^) ('f) = -A (k) if r = m 



= 0 if r + w, 

(11) 

and the expression (10) is then equal to 




(12) 

Thus the differential equations are satisfied 


Now considering the initial conditions, we put x- 

= 0, and find 


1 f ^ Pr8(y) 7 

A / \ 

27nJo A (k) 

(13) 

Also 

~ {^m (^) > 

(14) 

and 

2ti jo K A (k) ^ 

^nh 


2m /{ K ^ K A (k) 

dK (15) 


Now A (/c) IS a polynomial of the T^’th degree in /c, and JSrs (k) one 
of the (n - l)th degree When \k\is great enough the second term in the 
integrand is of order It therefore gives zero on integration, and we 
have 

^s = Us, (16) 

so that the initial conditions are satisfied 

Next, suppose that instead of the S ’ s being zero they are exponential 
terms of the form Since the solutions are additive it is enough 
to suppose that the initial values of the y’s are aU zero, and to consider 
the efiect of an exponential term in the first equation alone Then the 
equations are 


^ yi + ^2 ^2 + + 2/n ~ 

^ — /ji f 

( 5 + 1 ) 


(17) 



COMPLEX THEORY 


21 


The operational solution is 

A ^ (18) 

This IS to be interpreted as 

- Jl f {«) p 

(19) 

of “ the differential equations we find that the left sides 

of all vanish except the first, in consequence of the integrand Z! 
^aining as a factor a determinant with two rows identical! the first 

( 20 ) 

equations are therefore satisfied 

the inwln^ integrands in (19) are all 0 (<c-=) and 

S ncSr'? '^"’'"^ Thus all the ys vanish with i ^ 

thpn ^ foMtions that occur in physics can be expi essed by Fourier’s 

sJlntTrT® of exponentials, the proof that the 

olution can be carried out by operational methods is complete We 

shall see ater however, that Fourier’s theorems are not so co^enient 
to use as the formula 1 8 (9) based on the function H{a>) 

Tf wf ® integrals 2 1 (2), taken along the path L 

If we suppose the contour completed by a large semicircle from c + .oo 
to c - .00 by way of -oo , all the singularities of the integrand 7e Zh7n 

integral (1) aiound a large circle Now if a; is positive, and d, MU tends 

'' to 00, ?he integral 

iMvW becomes indefin- 

itely large, by Jordan’s lemma* Hence the integral alongX is equivalent 

to the integral around a large circle Thus if ^ is positm and 

^ W//c=0(/c~^) 

when K IS great and is positive, <l> Qp)H{x) is the same m 4>(p)l 

neltiv^ide^oflT’ semicircle on^the 

negative side of the imaginary axis is no longer an instance of Jordan’s 

e^rdT on Wer holds The integral around a large 

S o M °f tbe imaginary axis is, however, rl 

dumb e to the integral along Z, since by construction there is no 
singularity between these two paths If then </.(«)/« = 0 (/c-'O when k is 

Whittaker and Watson, Modern Analysis, 1915, 115 
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great, the limit of the integral around this large semicircle is sicto by 
Jordan’s lemma, and therefore <j} (j>) II (x) is zero 

Thus if 4> {p) IS expansible in descending powers of p, bognmuig with 
a constant or a negative power of^, we shall have 

= whena:>0'| qs 

^ (j>)j3'lx) = 0 when x co] 

The integrals 2 1 (2) for (x), where w is a positive integer, are 
divergent, hut these derivatives can be found by difleroiitiatimi Kvi- 
dently they are zero except when x = 0, when they do not exist 

Kgmx, dK = 11 {x + k) . 

This proves the result obtained by less satisfactory meaiiH in 1 8 
We can also translate into the form of a double integral the opera- ' 
tional expression for a general function From 1 8 (9) 

f{^)^r fQi)d{e'^^lI{x)], .. (:0 

J -00 

Inhere the integration* gives a function of p operating on II (x) Intor- 
pretmg by Bromwich’s rule 

/(^)=^ /I //(A) (0 

where the integration with regard to h is to be carried out first* 

If the function of p that arises when (3) is integrated has no Riiigu- 
lanty on the imaginary axis or on the positive side of it, and if ftirtlu^r it 
contains p as a factor, we can replace the integration with regard to 
in (4) by one along the imaginary axis Then we can put 

'< = ‘\ .. (5; » 

and we have 

/(j^)sink(x-A) dAdX .. ((!) 

This is Fourier’s integral theorem 

2.3. In practice the form 2 1 (2) is generally used m proforonce to 
2 1 (1) In problems involving a finite number of ordinary differoiituil 
equations of the first order the forms are equivalent for all poaitivo 
values of the independent variable, and since the independent variable 
IS usually the time, and we require to know how a system will behave if 
it IS in a given state at time 0 and is afterwards acted upon by known 

* The integral is of the type introduced by Stieltjes 

l' 
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disturbances, it la as a rule only poaitivo valuoa of tlio timo tluit con- 
coru iia 13ut when wo coiuo to deal with continiiona ayatoinH it usually 
happens that when tho oppiatimial solution ^ (p) m intoiprotod as ail 
iiitegial, the integrand has an infinit(> numlior of polos, and that no 
circle however largo can iiiclud(> all of them Hence tlie contour intogial 
2 1 (1) cannot be formed. But the hue inlogral 2 I (2) usually si, ill 
exists Again, it may happen that <l> (k) has a branch point at the origin 
or on the negative side of the imagiiiaiy avis Here again the conlour 
integral does not exist, but the line integral doo.s 

Ooiisoiiuontly in most of the wii tings of Iloavisido and Bromwioh, 
when a function of an operator oci'urs without the operand being 
stated explicitly, it is to be umlerstood that the opeiand is and 

that tho interpretation 2 I (2} is to be adopted ^riiis rulii will be 
followed in this work when wo come to treat continuou,s systeniH, and 
if it IS also supposed to bo adopted iii tho jirobloins of tho next chaiiter 
no harm will be done 

Considerable latitude is admissible in choosing tho iiath A Since by 
construction the integrand is regular at all points to tho right of A, we 
may replace A by any other path on its positive side ]novided it ('ndH 
at c + ico Again, the systems we aie to discush aio usually stahlo 
systems, that is, the jioloh of <f>(K) aie all on tho imaginary axis oi t/O 
tho left of it ^rium A can bo taken as near as we like to l,he imaginary 
axis A must not cross the imaginary axis it t.horo aro imles on tho 
latter, but it can ho taken to be a lino parallel to tho imaginaiy axis 
and distent i, fiom it 'I’lus is tho device most often used by Biomvvich; 
tho only advantege of tho piosont detmition ol A is that tho soluiioii 
is then adaptable to unstable systems as well as stable ones In all i-asos 
tho actual value of c does not atlect tho losults, so long as it is 
positive 


2.4. The principal operaf,ors involving branch-points ari' fractional 
powers of p Wo reiimie thou an intorpi elation o( p'^IlU), where « is 
fractional* Byouiiiilo 

/!"// U) . ..... ft) 

and the integial couvergi‘s if n I, A coiibrnr inehuliiig A ns jiart of 
itsoir, and such^ that the iiiti'graml is regular within it, is as shown, 
bvidoutly if .r is positive tho largo ((uadraiite make no oontrihiitnm. 

* Cf. Bmmwioli, Pruc Camh, Pint. ,SV . Ml (Uiai), 4U 127 
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The integral along L (= AB') is therefore in the limit equal and opposite 
to one along CBEF If further n is positiye the contribution from BE 
tends to zero with the radius of BE Thus A B, GB and EF contribute 
tbe whole of the integral Now on GB and EF 


K = ^ and /ig-" 

respectively, where ix is real and positive Hence on GB 
andon£!F ^ ' 


= - 6-!^ d/j. 


( 2 ) 

( 3 ) 



Therefore 

i e-i^ sin (« - 1) ^ 

__ sm nir r (n) 

TT ' ' 

But by a known theorem 

r («) r (l — w) = IT cosec MTT, 

and therefore 


p^Hisd) = 


nx—lh 

r(i^) 


( 4 ) 

( 5 ) 


( 6 ) 
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This result has been obtained for 0 < < 1 If we change n into — 
where - 1 < m < 0, we have 


r (m + 1) 

Powers of p outside this range can be found now by integration or 
dijBferentiation We have 






x^ 


( 8 ) 


r(^ + i) v(my 

which shows that (7) still holds when w — 1 is written for m (To justify 
this, the path AB must be replaced by FEDG ) 

Similarly 




rx 

y 0 - 


)V(m + l)^‘^ r(m + 2) 

By induction we may therefore generalize (7) to any value of m Hence 

= (10) 

Since when wi is a positive integer V{m + l) = m\ and when «2 is a 
negative integer P (m + 1) is infinite, the interpretations already adopted 
for integral powers of p are special cases of these 
In particular, since 

we have . 


P 


A. 


= 1 jt,f= ^ 

/J Tfx ’ 2 ij- 


and so on Also 


1 - 3^131 5 

-2 2 




( 12 ) 


(13) 


A lelated function that arises in problems of heat conduction is 
where a is a constant, and q denotes By Bromwich’s lule 

(14) 

On L the argument of is between ± Hence if a is positive the 
integral is convergent 

Immediate expansion of in a power series and integration term 
by term would not be legitimate, because all the resulting integrals 
after the first two would diverge But (14) is equivalent to an integral 
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along a path such as FEDG in Fig 1, and on this path we can proceed 
in this way Thus 

1 , a?K n , 

+ (15) 

All tlie positiTe integral powers of k give zero on integration The other 
terms are eq[uivalent to 


l—ap 


3 


3t 51 


Z\2\/xy ^ 2^5 \2jx^ 




’ZMX 

where, by definition, 

(11) 

By differentiation with regard to a we find 

( 18 ) 

We shall ^metimes need an asymptotic approximation to Erf w when 
w IS great We have 

1 - Erf w = —j-l dt = ^ f g-® u~idu 

\/Tr J w s/tt 

La 22 2221 


C~i)’^l 3 5 fco 

■ n/- ~T 2 2 2 L^”“-*<^+^>^.(19) 


on successive integrations by parts But the last integral 

rCO 

Jw 


2er^^ 


We have therefore the asymptotic approximation 


“ (l — Erf zt?) oj-L i_l 

VTT L 2 




2 2 2 2 2 


] ( 21 ) 
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CHAPTER III 

PHYSICAL APPLICATIONS ONE INDEPENDENT VARIABLE 

3.1. An electric circuit contains a cell, a condenser, and a coil with 
self-induction and resistance Initially the circuit is open It is suddenly 
completed , find how the charge on the plates varies with the time 
Let y be the charge on the plates, t the time, K the capacity of the 
condenser, L the self-induction, and R the resistance of the circuit, 
and let E be the electromotive force of the cell Write <r for djdt 
The current in the circuit is y, and the charging of the plates of the 
condenser produces a potential difference yjK tending to oppose the 
original e m p Then y satisfies the differential equation 

E-^=Ly + By ( 1 ) 

Initially y and y, the current, are zero Hence the subsidiary equation 
is simply 

(La^ + B(T + -^y = E, (2) 

and the solution is 

E 

T 

L(t^ -h Rcr -h — 

j6l 

If now the denominator is expressed in the form X (o- + a) (<r + yS), the 
interpretation is, by 1 5 (9) 

Since a + /3 and a.p are both positive, a and p must either he both real 
and positive, or conjugate imaginaries with positive real parts In either 
case y tends to KE as a limit, as we should expect 
We notice incidentally that if the circuit contained no capacity or 
self-induction the differential equation would be simply 

B<Ty = E ( 5 ) 

Hence if a problem has been solved for simple resistances, self-induction 
and capacity can be allowed for by writing io- + i? + -A for ^ For this 

Ecr 

reason this expression is sometimes called a ‘resistance operator,’ and 
the method generally the ‘method of resistance operators ’ 
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physical applications 

3.2. The Wheatstone Bridge method of determining Belf-Inducti on, 
^ this method the unknown inductance is placed in the first arm of a 
Wheatstone bridge, the fouith arm is shunted, a known capacity boinc' 
placed in the shunt* 

First consider the ordinary Wheatstone bridge, the rosistaucos of the 
arms being 5i, R^, R^, R^^ let ai be the current in Mi,y that in /4, ^that 
through the galvanometer, and b the resistance of the battery and loads 



Iko!- Biy+Gg = 0, 

R»^-Biy-{Ri + R^+G)g = 0, 

^i«> + y) + Biy + R4.{y+g) = E, 
and on solving we find 


.. ( 1 ) 
• ’W 

..(3) 


x + y 

w + ±i, + R, + R,) + + R^) (*) 


R^ Rs — Ri R^ 

' * — ‘ ■ -"V -r 

it b IS large compared with R, and R^, we have nearly 

a!+y = Elb, 

and = (.EIb){R,R,-RE.'\ 

® ^^1 + -ti, + M, + E,) + {Rf+R^) 4.7^^) ■ 


( 5 ) 

( 6 ) 


trfitTir^ W n ""r self-induction m 

WW T ^ ^ arrangement in the 

fourth arm be as shown (Pig 3) The resistance of the mam wire is Z 



IPjg 3 

* Bromwioh, PUl Mag 37 (1919), 407-419 


Further leferenoes are givea 
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that of the ahiuited portion of it r. Huppogo the Hhnnt to have a 
resistance iS' Then the effective roaistanco of the whole arm would bo 

It, ■ ■ (7) 

If the shunt contiuim a capacity h' wo muat replace xS'by A'+l/AV 
Hence in the formula for g m luiiHt replace /f, hy Ao- and by 


//. 


» A'<r 


(p I iS')AV+ I 

The result expressoH the current through the galvanoinotoi when the 
battery circuit is suddenly closed. 

It can be shown that in actual conditions q cannot vanish for all 
values of the time A sufficient condition tor this would bo that the 
operator should bo identically zero, then g would bo 

identically zero whatever the remaining factor might lopiosont But 
with our modifications this factor becomes 

(«) 

Multiplying up and equating coefficients of powers of cr to zero, wo 
find ’ 

/ACr I *0 - (<J) 

— LRi f (/4 /fj //, //<)(? +aS') K I- Rit^K .(), {^10^ 

R'iRn — RiR^ - 0 . . , _ (11^ 

From the construction of tho apparatus r< //„ r<.r + H a’hus (9) (mui 
hold only if r Mi and iS' (). The shunt wire must bo attached to the 
ends of R^ and must have zero resistance flubstitufang in (to) wo find 

L R^Rih', .. .. • (12) 

together with tho usual condition (11) for permanent balance of tho 
bridge 

Actually these conditions for (omplete balance cannot bo completely 
satisfied, but for the dotoimmation of A it is not necessaiy that they 
should fciuppose the galvanometer is a hallisfuc one, and that if is so 
adjusted that there aio tie jiermanont current and no throw on closim/ 
the circuit Thus ^ 

LllUf/ 0, f gtU 0 . 

/ jr» J ^ 

If 17 18 exjiresHod in the foim/((r)/A'(tr), wo have 

/(<») , ./(a) 


(liO 


it 
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where all the a’s, lu the conditions of the prohloin, will havo lu'gativo 
real parts Then the condition that g shall tend to zero f,nvos 

/(()) = () .. (lAj 

Also gdt = % •(<«) 

Equation (15) shows that the operational form of g contains tr iw a 
factor Also we can write 




/(“) 


(17) 


aF'ik)cr-a 

SO that (16) IS the limit of (— gjcr) when a- tends to zero. The vaiusliinj^ 

r<x> 

of Lim g and gdt imply therefore that the operational form of g 


contains as a factor Hence the modified form of /4 7/4 iniwi 

contain cr^ as a factor Then (10) and (11) still hold, (0) no longw 
holds We now have 


R^Hi^ — Oj . , . ( I H) 

LBi^lt^r^K, .. (I'J) 

which gives the required rule for finding L 


3.3. The Seismograph In pnnciple most seismographs are Euler 
pendulums— pendulums with supports rigidly attached to the earth, so 
that when the earth’s surface moves it displaces the point of support 
horizontally and disturbs the pendulum The seismograph differs from 
the Euler pendulum as considered in text-books of dynamics in two 
ways instead of being free to vibrate in a vertical plane, it is con- 
strained to swing, like a gate, about an axis nearly, but not quite, 
vertical, so that the period is much lengthened, and fluid viscosity or 
electromagnetic damping is introduced to give a frictional term pro- 
portional to the velocity The displacement of the mass with regard to 
the earth then satisfies an equation of the form 


x+2Kx + r^a, = \^^ 

where i is the displacement of the ground, and s, and X are constants 
of the mstoent Suppose first that the ground suddenly aequinis 
a fimte velocity, say unity Then 



* Some instruments, such as that of Wieohert, are not on the pnnoit.le of the 
Euler pendulum, but nevertheless give an equation of this form. 
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and our aubsidiary O(j,uation is 



(er® 4 - 2 /«t + n^) x =5 ktrll (f) 

(3) 

Put 

0*^ + 2^0- + 71? - (o- + a) (a + p) 

(4) 

Then 

(c + a) (o* + p) 

(5) 


= 0 when t<Q . 

(6) 

and 

=• when t>0 . 

(7) 


Tho recorded displacemout x theroforo begins by increasing at a finite 

1 


rate A, roaches a maximum X ** ^ softer a time 




and 


then tends asymptotically to zero 

If a and ^ are real, and /3 less than a, we see that the behaviour after 
a long tune depends mainly on now as the experimental ideal is 
to confine the effects of a disturbance to as short an inteival after wax ds 
as possible, we see that wo should make p as large as possible But 


yS- — 

K 1 


( 8 ) 


and for a given n, p is greatest when k - n This is tho condition for 
what IS called ajieuodicity , the roots of tho period cciuation are equal, 
real, and negative Many seismographs are at ranged so as to satisfy this 
condition Tho soluiaon is then 


"^0 when t 'U) (10) 

-kU when i>i) (n) 

# 

Tho maximum displacement is now at time l/n after tho start, and is 
ecpial to kjen 
If K<n, we can put 

. ( 12 ) 

Then tho solution is 

7 


( 13 ) 
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and the motion does not die down so rapidly as for k n. The aiionodic 
state therefore gives the least motion after a long interval for a given 
value of n 

In practice, however, k is usually made rather loss than n ; in the 
Milne-Shaw machine, for instance, k is about 0 Tit ^I'he motion is then 
oscillatory, but the ratio of the first swing to tho second is or about 
20 But X vanishes after an interval r/y from tho start, or about 4//i, and 
ever afterwards is a small fraction of its iirst maximum 'riu* reduced 
damping effect after a long time is considoied loss important than tho 
quick and complete recovery after tho first maximum 'I'ho time of tho 
first maximum is 1 l/n from the start, as against l/« for tho aporiodio 
instrument and 1 57/ra for the undamped one 



3.31. The Galitzin seismograph is similarly arranged, hut tho motion 
of the pendulum generates by electromagnetic induction a current, 
which passes through a galvanometer If x is the displacement of tho 

pendulum, and y that of the galvanometer, the differential cmtatioiiH 
are 

«+ 2Kia! + %“<» = Xf, ,, 

y + 2Kay + maV‘=/^ (2) 

Supposmg both pendulum and galvanometer to start from rest, wo have 



(o^-i-2Ki<r+rei‘')(o^ + 2Kjw+mj“j ' ‘ ‘ 'w 
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As a rule the two interacting systems are so arranged that k and n are 
the same for both, and both are aperiodic, so that K-n Then 




(o- H- ny 


(4) 



If we suppose as before that 

(5) 

we have y = (t), 





(6) 


The indicator therefore begins to move with a finite acceleration, instead 
of with a finite velocity as for the pendulum The maximum displace- 
ment follows after time (3-^/3)/y^=l 27/w, the mirror passes through 
the equilibrium position after time 3/w, and there is a maximum dis- 
placement in the opposite direction after time 4 73/w The mirror then 
returns asymptotically to the position of equilibiium The ratio of the 

two maxima is / + \/3)^ = 2 3 In comparison with an instrument 

such as the Milne-Shaw, recording the displacement of the pendulum 
directly, the Galitzin machine gives the first maximum a little later, 
the first zero a little earlier, and the next maximum displacement is 
larger in comparison with the first maximum 


J 


3 
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In an actual earthquake the velocity of the ground ib annulled by 
other waves arriving later, the complete motion of the sewmograph is 
a combination of those given by the sepaiate displaiiomontH of the ground. 

3.4. JSesmmce A simple pendulum, originally hanging m ttiinhhi mm, 
IS disturbed for a finite time by a force varying harmonically in a period 
equal to the free period of the pendulum Find the motion afUu the 
force is removed 

The differential equation is 

a! + n^a!=/mint, . (j) 


-/ 


a; =f j 


0-* + II? 
ficr 


The solution is then , 

((r“ + «“)■•*’ 

noting having to be added on account of the initial oouditioiiH 
To evaluate this, we notice that 

Differentiating with regard to n, we have 

2wcr _ 1 

(<? + f?)^ 




and 


f 

® ~ 2 ?? ~ <^ 0 ® ^0 


(-0 

(») 


(-0 


(5) 

(«) 


Suppose the disturbance acts for time ^Trjn, where r is an intogor. At 
the end of this time 


(7) 


a;=() 

The subsequent motion is therefore given by 

^ ^-^^^COS(w!f-nr), 

rirf 

• .. (8) 

to!th^^!£ffr aroattodied 

particles of mass!! fixed One of the 

LquentmotiTofrt ? transverse impulse J Pmd the sub* 

Ifr r* Lp rt!7 f ^ ® (Ii^tercoUegiate Examination, 1923.) 

and nC'dXc! J" tension, 

..y th. ,h«e ■" 
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Xi — A, (^ 2 xi 

( 1 ) 

1^^272 = - X + 2X2-X3), 

(2) 

~ X ^ OD2 2^3)5 

( 3 ) 

P 

where X = — ^ 

ml 

( 4 ) 

Initially all the displacements are zero, X2 = X3=^ 0 , 

Hence the 


subsidiary equations are 



(or^ + 2X) Xx - ^^^2 = o-Ijm, 

. (5) 


— ^Xi + + 2X) X 2 ” Xi??3 = 0, 

(6) 


~ \X 2 + (cr^ + 2X) X 3 -O 

- (7) 

Hence 


(8) 

and 

/21 , 2X2 V x<7 I 

<r2 + 2xj^=‘^<72 + 2Xira 

(9) 

Multiplying by + 2X we find 


and 

(7<t2 + 4X) (3tr2 + lOX) ^2 = 20\o-Ilm, 

(10) 

20X<r J 

'®'"(7a^ + 4X)(3<T2 + 10X)OT 



_20I r 7cr So- ) 



58 170-2 + 4X 3(;^+ioXj 



_10/ fsma# 

29wt a P 

(11) 

where 

a2=*X, /32 = ^X 

(12) 


We notice that the mode of speed s/(2X) is excited, but does not need 
to be evaluated because it does not affect the middle particle 

3. 6- Hadioactzve Disintegration of TD amum The uranium family of 
elements are such that an atom of any one of them, except the last, is 
capable of breaking up into an atom of the next and an atom of helium* 
The helium atom undergoes no further change The number of atoms 
of any element in a given specimen that break up in a short interval of 
time IS proportional to the time interval and to the number of atoms of 
that element present If then w, x^ are the numbers of atoms of 

the various elements present at time t, they will satisfy the differential 
equations 

* We neglect jS-ray products, for reasons that will appear later 

3“2 
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du 

dt 




\ 


dSTi 

\u 




d^% 

dt 




( 1 ) 


dx„ 

dt 


n, l 


Suppose that initially onlyiiramum ih promuit Tlitin wlnni t 4), w ?%, 
and all the other dependent variahloH are zeio. 1'hen the Hiihhidiary 
equations are 

((T-f k)u <r/^at 
(cr + Kj) /T'l 
(cr-fKa)^a 


(o* + -I ^n, SI *^V t» 

o*/r ^ K- 


The operational solutionB are 


(TU^ KtTHa 

0*+k’ (<T k) (or I* ^ 


KKi (THo 

((r 4* k) (<r K,) (cr B Kjj) ’ 


.( 1 ) 


(0" + K)(0r4Ki), j) J 

These are directly adapted for interpretation hy the part4al-fmt*iioti 
rule In fact 




Ki ■**“ K 


^2 — KKi^q 


■ 1 

■(ki-k)( 


(kj - k) (kj - (f) "'’ (k - 'C,) (/c/. 


fi Xlf 


(*-K3) 


ft ; 


} ...(4) 


Xn = Uo~ 
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Of all the decay constants k is much the smallest If the time elapsed is 
long enough for all the exponential functions except to have become 
insignificant, these results reduce approximately to 

= -^-*0 ( 5 ) 

Ki ^2 

With the exception of the last, the quantities of the various elements 
decrease, retaining constant ratios to one another 

On the other hand, if the time elapsed is so short that unity is still 
a first approximation to all the exponential functions, we can proceed 
by expanding the operators in descending powers of o- and interpreting 
term by term Hence we see that at first will increase in proportion to 

t, Xi to f, and x^ to T /.a 

In experimental work an intermediate condition often occurs borne 
of the exponentials may become insignificant in the time occupied by 
an experiment, while others are stiU nearly unity We have 


ODf' 


or 4- /Cy 

and if K,t IS smaU we can neglect the second and later terms in com- 
parison with the first Thus in this case 

-1 

If Xr.-y IS of the form f, we can put 


Xr-- 






Q-Kft 


(«) 


■ cr*'^^ or + Kj. or*"^^ 

If K,# IS small, wo can replace the exponential by unity and confirm 
(7) If It IS great, 

= I e~‘^rtdt = - + 0 (e-*’-'') 


and on continuing the integrations ^ ^ 

= or-* — = — 7; 


Thus 


Xr = 




( 9 ) 


' -y- # 

Classifying elements into long-lived and short-lived according as «,.# is 
amflll or large for them, we find that the quantity of the first long-lived 
element after uranium is proportional to t, the next to f, and so on 
All |8-ray products are short-lived when t has ordinal y values 
Eadium is the third degeneration product of uranium In jock 
specimens the time elapsed since formation is usually such that the 
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relations (5) have become established As a matter of observation the 
numbers of atoms of radium and uranium are found to be in the con- 
stant ratio 3 58 X IC"’ This determines k/ K g Also the rate of break-up 
of radium is known directly in fact 


1/ Kg- 2280 years 

Hence 1/k = 637 x 10® years 

This gives the rate of disintegration of uranium itself 
A number of specimens of uranium compounds were carefully freed 
from radium by Soddy, and then kept for ten years It was found that 
new radium was formed , the amount found varied as the square of the 
time This would suggest that of the two elements between uranium 
and radium in the series one was long-lived and the other short-lived 
Actually, however, it is known independently that both are long-lived 
The first, however, is chemically inseparable from ordinary uranium, 
and therefore was present in the original specimens , initially, instead 
of a!i = 0, we have 


K 


0^1 — 'Wq 


For the next element, ionium, we have 


the variation of oox being inappreciable in the time involved Also 

Xz = K2Cr''^CCz = ^KK^U^f 

Soddy* found that 3 kilograms of uranium in 1015 years gave 

202 X 10-“ gm of radium Hence, allowing for the difference of atomic 
weights, 

W«o=7 1xlO-“ 

'‘a = 8 64 X 10-7y®0J^, 1 /^ 2 = 1 16 X 10® years 

This gives the rate of degeneration of ionium Soddy gets a slightly 
lower value for l/xa from more numerous data 


3.7. Some dynamical applications Suppose we have a dynamical 
system specified by equations such as those in 1 72, save that sc is 
replaced by t, and that the system is at first at rest and then disturbed 

by a force 8r applied to the coordinate y. The subsidiary equations 
are 

= 0 (m=t=r), \erM'y»-8r {m^r) (1) 

* Phil Mag (6) 38, 1919, 483-488 
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Writing A for the determinant formed by the e’s, and for the minor 
of Brs in this determinant, we have the operational solution 


2 /.= 



( 2 ) 


If the determinant A is symmetrical, so that 

Ers = Esr (3) 

we see that a given force Sr applied to the coordinate yr will produce 
precisely the same variation in as the same force would produce in 
y, li it was applied to ?/, Thus we have a reciprocity theorem applic- 
able to all non-gyroscopic systems, linear damping does not invahdate 
the argument if the terms introduced by friction contribute only to 
the leading diagonal of A 

If the forces reduce to an impulse, so that Sr can be replaced by o-/,., 
the solution becomes 

= (4) 


Wo can evaluate the initial velocities by expanding in descending 
powers of o- The first term is 

y,^^\jr ( 5 ) 


where A is the determinant formed by the a’s, and Ars the minor of 
a„ m it Hence the initial velocities, found by operating on this with 

or, are 

<ry. = ^Jr • ( 6 ) 


Now the constants are merely twice the coefficients m the kinetic 
energy, which is a quadratic form Hence the determinant A is sym- 
metrical whether the system is gyroscopic or not, and the reciprocity 
theorem for impulses and the velocities produced by them is proved 
The subsequent motion can be investigated by interpreting accorffing 
to the partial-fraction rule But let us consider the simple case where 
the system is non-gyroscopic and frictionless Then 

A = JLn (cr“ + a!‘), • ('i) 

where the a’s are the speeds of the normal modes Ihen 




sin oit, 


... ( 8 ) 


where Er,(ro?) denotes the result of putting - a’* for cr“ in Er, The 
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contribution of the a mode to the initial rate of change of y* is therefor© 

„ _ (- a“) J. 

“ (dA/da^),.^y^ •• ( 9 ) 

An immediate consequence of the presence of the factor iSv„ (- a“) m 
the numerator is that if an impulse is applied at a node of any normal 
mode, that mode will be absent from the motion generated 
Another illustration is provided by Lamb’s discussion* of the waves 
generated in a semi-infinite homogeneous elastic solid by an internal 
disturbance The normal modes of such a system include a typo of 
waves known as Eayleigh waves These may be of any length, and 
involve both compressional and distortional movement, if the depth is 
the amplitude of the compressional movement in a given wave is 
proportional to and that of the distortional inovonient to 
where a and /3 depend only on the wave-length Lamb found that if 
tbe original disturbance was an expansive one at a depth/, the ampli- 
tude of the motion at the surface contained a factoi but if the 
original disturbance was purely distortional, the corresponding ampli- 
tude contained a factor e-^ These factors are the same as would occur 
in the compression and distortion respectively at depth /in a Eayleigh 
wave with given amplitude at the surface 


OHAPTEE IV 

WAVE MOTION IN ONE DIMENSION 

4.1. In a large class of physical problems wo meet with the difleren- 
tial equation 

Sf dsc^~^’ — ( 1 ) 

where # is the time, so the distance from a fixed point or a fixed plane, 
y the independent variable, and c a known velocity Let us consider the 
solution of this equation first with regard to the transverse vibrations 
ot a stretched string In this case we know that 

where P is the tension and m the mass per unit length Wnte o- for 
S/3#, and p for djdso Suppose that at time zero 


y 


-/(-), |: 


Phil Tram A, 208, 1-42, 1904 


( 3 ) 
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where / and F are known functions of x , that is, we are given the 
initial displacement and velocity of the string at all points of its length 
Then we are led by our previous rules to consider the subsidiary 


equation 

(cr‘^ - c^p^) y = a^f (x) + cr P (x) 

(4) 

or 


(5) 


or*^ 1 / ^ ^ \ 


But 

a^ — C^p^ ^\cr—cp cr + cpJ 




(6) 


or-* — 0^ \(r-'Cp or + cp) 

2c^ ^ p 

In (7) the Xfp has been put last in consequence of our rule that opera- 
tions involving negative powers of p must be carried out before those 
involving positive powers Hence 

= 4 {fip + +/ 


_f_i?’(a,) = i ^ P{x)dx 

(T^-C^P ^ ' 26 Jo 

r— rflo + ct p X ct T 

-iU. 

. (9) 

Sic j IS — et 


and tlierefore 

2/ = i {/(« + 0*) +/(® " 



( 10 ) 


This IS D’Alembert’s well-known solution 

This cannot, however, be the complete solution Equation (1) Has 
been assumed to hold at all points of the stnng , but if any ex ema 
forces act these must bo included on the nght of the equa,tion of motion 
In such problems the ends of the string are usually fked, and ructions 
at the ends are required to maintain this state, in the complete equa- 
tions of motion these reactions should appear They are unknown 
functions of t, and therefore necessitate a change in the mode ot 

solution 
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But D^emberts solutloa can be adapted so as to fulfil all the con- 
ditions We notice that the solution consists of two waves travolliiiff in 
opposite directions with velocity c If the initial disturbance is confined 
to a region within the stnng separated by finite intervals from both 
ends It will take a finite time before D’Alembert’s solution gives a 

bSar^* reciuired to maintain tho 

the ™ ’ and the solution will hold accurately until one of 

me waves reaches one end 

lenvfh nf disturbance is specified only foi points witliin the 

by the last paragraph D’Alembert’s solution would bo couipleto if the 

T- stretching 

from oo to + CO , and disturbed initially so that/ (<r) and FU ) are both 

anhsymmetrical about both .= 0 and ir ^ f, y will Lsh evoViftei It 

fZJS Ti T “*<■■>. “'i 

q ivalent to that for tho actual string from a; = 0 to a; - ^ Thus wo have 

veTocitvS d the actual values of tho initial 

velocity and displacement between ai=0 and but with the dis- 
placement and velocity outside this stretch so spocifiod as to bo 

SfonTr^re\rf ^ 0 and a: = f , Li D’Alombort’s 


_ (T^ 


^ O'® 


•• (I) 

except that the effects of the initial velocity will not bo considered at 

Scion * as If a were a constant 'tIi! 

j stifieation of this procedure is given later in 4 6 The boundarv con- 

variation of parameters, we assume that the solution is 


y=A cosh ~ +jB sinh — , 
^ c 


where A and £ are functions of a; subject to 


A'cosh^ + ^smh‘^'® = 0 

6 0 


( 2 ) 


( 3 ) 
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Substituting in the equation (1) we find 

A' sinh — + cosh — / (x) (4) 

c c c 

Hence 

A' = ^/ (^) sinh ^ , F = - {p) cosh™ (5) 


Now 3 / must vanish for all time when , hence A (0) — 0 Thus 



A=r^/(i)swh^di 

(6) 

Also 3 / must vanish when x=l Hence 



B (1) smh j = -A{t) cosh 

(7) 

and 

B(t) = - coth - f -/ (1) smh - dl 

(8) 

giving 

B(x) = - coth— [ —/(i) smh y di 

C Jo 0 0 



+ £‘^/(i)oosh^di 

(9) 

In all we find 



^=j 

j ^ / (^) smh ^ ^cosh ^ - coth — smh — ^ di 


+ 

j 

[1/(0 smh =(eosh ^-c»th f s.nh f) di 





+ 


(10) 


Before carrying out the integration with regard to we can interpret 
the integrand by the partial-fraction rule Each integral vanishes when 
or = 0, and therefore no constant term appears in the solution But 


sinh dljc vanishes when 

( 11 ) 


where r is any integer, positive or negative Hence 

f- 


o* - 0-^sinh (T (^— i2?)/c I I 

- smh ^ 


criTT 


7 


cosh (tltt) 


I 


0ri.irGtll 
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(,3) 

This IS the solution given by the method of normal coordinates Puttinir 
t zero, we obtain the Fourier smo series ^ 

y(^) = 7p(^)su/f- 


?7r.r 

Hill ~^~dc 


The effect of an initial velocity can bo obtained bv 

OfifiSS */ 


process 


(H) 

a Hiinilar 

Bat m pnwtiee the sohtion by ti iBOimuiotrical mucs a not often tlio 
teott convenient fom. It nennlly oonv,i„«, „b„|y, J^Xt T mote 
eenone that .te term ansgeet, little about tlio Mtiito Inlroitml 

oetoial form of the string at any iiiatont it le tiocosaary to find eoino wav 
of summing the senes, which may bo rather dilheiilt i inorrc3iteS 
method IS often the following > To have . 00 , . that tlm 1 Z “^3 of 
any operator vijid for positive values of tho argninont ", wZZ te 

i -6” 


) (If') 


8inhcr//c * i +«-2<rl/0+ ^ 

and when interpreted by Bromwich’s rule this will give nee to a HeWnf, 
m powers of e-« the integrand Bnt smeo the roa Z m 

ZZZn™ I" “ “tsolntoly convZni 

Zd IZZtoZ this Z” “““ ■* “ 

tee operator in this way and to interpret term by term Then for 

sinho-^/cfiinho-ff-a,)/o 

siiho-z/ir ~ (1 - 

andfor^x., d " 

sinh crmlo sinh <r (I- ^)fc 

ie-''(f-»)/r (1 _ s-i.r^K') 

+ N /,„a, 

exponentmkoTt^^^^^^^^^^ 

* Heaviside, MlectroviagmUo Theory, 2, 108-114 
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Then ^/(i) e-<^ di = J/(i) ( 18 ) 

But according to our rule is zero wheni^-(a-i)/c is negative 

and unity when this quantity is positive Hence the integral (18) is 
zero unless lies between and 4, and then it is equal to / («^ - ot) 

Similarly (19) 

and e-*j-(^-a)/e = 0 li ct + a 

=- 1 if ^ < cif + a, 

BO that ( 19 ) IS zero unless ct^-a lies between 4 and 4, and then is 
equal to/(ci^ + a) 

It follows that at any instant the effect of the initial displacement 
at the point ^ is zero except at a special set of points where one 
or other of the quantities a ± cif is equal to i Since a does not involve 
4 we see that this way of expressing the solution reduces it to a set of 
•waves moving in each direction with velocity c The first three factors 
in (16) give 

1 {x-§)lc ^ {x’^$)/c — ^-o-{(a5-f)+2 {l-x)}ic ^ g~(r{(a?+f)+2 (Z-a7)}/cj. 

( 20 ) 

for values of x greater than i The first term gives at time 

(a}-i)jc, the second -i/(l) at time ix + $)lc, the third -|/(^) at 
time {x-i + 2 (l-x)}/c and the fourth ^/(l) at time {x + i + 2 Q- x)]l c 
The first term represents the direct wave from ^ to x, the second the 
wave reflected at ir = 0 , the third that reflected atx = l, and the fourth 
one reflected first at a; = 0 and then again at x= I The term in 

the last factor of ( 16 ) will give four further and similar pulses, each 
later by 2 llc than the corresponding pulse just found These are pulses 
that have travelled twice more along the string, having been reflected 
once moie at each end Similarly we can proceed to the mterpretation 
of later terms as pulses that have undergone stiU more reflexions 
The part of the solution arising from values of ^ greater than x may 
bo treated similarly The interpretation in terms of waves is exactly 
the same 


4.3. As a particular case let us consider a string of length I, origin- 
ally drawn aside a distance 17 at the point x = b, so that initially it lies 
in two straight pieces, and then released Then 

ya-'nxib 0 <x<b\ 

ya = fj(J—x')j(l—b') b<x<l) 
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and the subsidiary equation is 

0-“ 

die e 

The solution that vanishes at the ends is 

y~tl^ + A siiih sinli h) 


0 .r < h 


y = A aiiih Hnli b^i<l 


( 2 ) 


C3) 


where the constants have boon chosen so as to niako y continuous at 
a>-l Also a discontinuity in dyjdoo at this point would unjily mi infinite 

FT'* 

dy/dx snail be continuous, which gives 

S + JTi) + {cosh j sinh f>) + sinh cosh J (-?-/>)} 0 


and on simplifying 


^ ^ t O’/ 


•• ( 4 ) 
■ (5) 


.( 6 ) 


Thus 

b(l-b)\ I or sinh(r//c 

b (I b) u o- smh crtjc / ^ ^ J 

Interpreting by the partial-fraction rule, we notice that tho contnbu- 

nidepondent of <r, while tho rest 

gives 

^ b(l---b) r^i I ~ OOB 0<X<1 (7) 

^Wtody. w, caa mtopyet tho .otaon for 0 < * < J .„ oxpoMontml,. 

Sinh (Tx/c smh cr (I In 

= |c-o-(S-»)/r- (] « 


and 

Then 


(l 6 (l + 4 . jij’*' W/f* ^ ^ ^ 


- ss C/ Jf > 0 

<T * ** 

6--i^-^)ht = ot~{b-x), ... 

Ea7leigh, Theory of Sound, I, 1894, 185 


•• (9) 

-(H)) 
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when this is positive, and otherwise is zero Thus y retains its initial 
value yjaill unaltered until time (6 - «)/c, when the wave from the part 
where » > 6 begins to arrive After that we shall have 

f-X i ~ ^ = KCT) { G ^ 

(11) 

This solution remains correct until ct = b + a}, when the wave reflected 
a.tx = 0 begins to arrive Thenceforward 

— e-o'(6+»)/«) — = {ct — (b— ic)} — {ci — (b + x)} = 2x (12) 




Fig 6 
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TLe part of the string reached by the first reflected wave is therefore 
parallel to the original position of the part where h<x<l 

When + h), a wave reflected = I arrives, and 

+ + (14) 

This holds until ct = (b 




(15) 


— in the next phase 
/- b\ y)CC 

'b(l-bf[^~l) = T 

When the whole of the string is back in its original position, 

the term m then begins to affect the motion, and the whole 

process repeats itself We see that at any instant the string is in three 
stoaight pieces The two end pieces are parallel to the two portions of 
string in its initial position, and are at rest For the middle portion 
e gradient dyjdx is the mean of those foi the end portions, and the 

transverse velocity is + 2 j The middle portion is always either 

extending or withdrawing at each end with velocity c [Fig 6 ] 

4.4. A uniform heavy string of length 2/ is fixed at the ends A 
particle of mass m is attached to the middle of the string Initially the 
strmg is straight An impulse J is given to the particle Find the sub- 
se^quent motion of the particle (Rayleigh, Theory of Sound, 1, 1894, 

The differential equation of the motion of the string is 




= 0 


( 1 ) 


ont'rt symmetry we need consider 

0 < * < ^ Suppose that the displacement of the 
« ^ ^ When t-Q, y and are zero except at «: = 0 The 

isIherSe " ^ The solution 

sinW(£-^c 

smh aljc (^) 

" ‘‘‘ of ““‘“I »f *I>9 


m 




2P 


\W, 




If P be the mass of the string per unit length, 

P = pc^ 




( 4 ) 



WAVE MOTION IN ONE DIMENSION 49 


On account of the initial conditions the subsidiary equation for the 
particle is 



= - 2P77 - coth — + Jot 
c c 

(5) 

Hence 

^ m<TC + 2P coth crljc 

(6) 

Put 

^ ^ _ i^siss of half the string 
m mass of the particle 

(7) 

Then 

1 

II 

(8) 

and 

Ujme 

^ {o-ljc) + h coth (cr//c) 

(9) 


To interpret this solution by the partial-fraction rule, we recollect that 
the system is a stable one without dissipation, and therefore all the 
zeros of the denominator are purely imaginary Putting 


crljc = t(0 (10) 

we see that the zeros are the roots of 


co = ^COta) (11) 

There is a root between every two consecutive multiples of tt, positive 
or negative, and the roots occur in pairs, each pair being numerically 
equal but opposite in sign Then 


IJ ^ I 

yj = ^ ( ^ 

. me LUic (1 + k cosec® w) (l/c) 


yKOOtll 


( 12 ) 

the summation extending over all positive and negative values of o>, 


2iJ 

= — S ^ gin (13) 

me o) (1 + ^ cosec® co) ^ ^ ' 

the summation being now over positive values of w 

If a root of (11) IS w-TT + A, where n is great and A < tt, we have 

(nw + A) tan A = A (14) 

and X=klmr (15) 

approximately The series (13) converges like :§ ~ sin® co, or 5 A It 

(0 W 

therefore converges fairly rapidly at the beginning, but more slowly 
later Also the w's are incommensurable, and therefore the motion does 
not repeat itself after a finite time, thus the labour of computation 
would be great 


j 


4 



50 


WAVE MOTION IN ONE DIMENHION 


To express the solution in ttmns of waves, it is oonvoniont to elmiige 
the unit of time to //o, tho time a wave takes to travt'rho half the length 
of the string Also J/m may bo repliwetl by V 'riien 

V I'd e *") 

~ o- + coth <r (<r 1 X) (ir 


Cl'- 

r /, -a 2K- o„„ 

<r + A \ <r I / («■ 1 kj' ) 


do) 


(ir 1 k)' 

The first tonn is zero for negative values of the tiino, for jKtsitive 
values it IS eipial to 

[(1-0 '0 (n) 

After time 2 tho second tonn no longer vaiiisluss We have 

k 1 

(o- + /)“ k kitr t i) (ir I / 1** 


ti< ** 


.(18) 


1 i 

k a ' 

and the second term is oiiuivalont hi 

- {1 - e - A (if 2) e '■*»! for t. 2 (I a) 

The third term is zero for ^ < 4 ; for gi eater valiu's it is easily found to 
be 


2F 


+A(< 4)4 r(« 4)ne*i**>l- 


.( 20 ) 


The process may be extended to detorniine the motion up to any in- 
stant desired The entry of a now term into tho solution corrosponds 
to the arrival of a new wave reflected at tho onds. 


4.6. A uniform heavy bar is hanging vertically from ono end, and a 
mass m is suddenly attached at tho lower end. Find how tho tension at 
the upper end varies with the tune. (Love, h'lemtkitif, S 28.T.) 

If w be the distance from the upper end, and t/ tho longitudinal dis- 
placement, 1 / satisfies the differential wpiation 

A- 
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where p is the density, JS Young's modulus, and JP the external force 
per unit volume, in this case pg Put E/p = c^ and let the displacement 
of a particle before the weight is attached he y,. Then 

( 2 ) 

When = 0, yo = 0 ; and if the bar be of length I, = 0 when ^ 


Hence 




After the weight is attached, we still have 

dx^ ^ da^ ’ 

while when t = 0, y = y^ and 9^/3# = 0 Hence the subsidiary equation is 


<r^y — c 






and the solution that vanishes with x is given by 

y-y^^ A sinh <txIc, (6) 

where A is independent of x 

If ‘sr be the cross section of the bar, the equation of motion of the 


mass m is 


d^y dy 

= (7: 

the derivatives being evaluated at /r = Z The subsidiary equation is 

mcr‘‘y=mg + mcr^y,-ma^, (g; 

which gives on substituting for y from (6) 

(<j^ sinh — + cosh — ) 4 ^ Zg] 

The tension at the upper end is 

gpl + JEcrA jc = gpl + q 

Jbj'm cosh {(rljc) + moar smh {(rile) 

If h be the ratio of the mass of the weight to that of the bar, 

j m me A , ^ 

P'csV Em" c' 


4-2 
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and the tension is 


gm 


( 12 ) 


\_k COsll((r//c) + k (cr^/0 Smh (rl/cj 
We see that gmj'mJo is the tension due to the weight of the bar alone, 
and grrijm is the statical tension due to the added load alone To evaluate 
the actual tension, we expand the operator in (5) in powers of 
Taking Ijc for the new unit of time, we have 


]c(T sinh <T + cosh o- “ (Jc<x + 1) - (k(T 
26^°^ r io 

+ 1 L ^0" + 1 


■ 1 ). 

-1 


y- 2 (r 


j-2<r 4. 


/Ao-- iY_ 

W+iy ” 


■4cr 4 , 


(13) 

The first term vanishes up to time unity, and afterwards is equal to 

2 (1 - e-«-i)/*) (14) 

This increases steadily up to time 3, when the next term enters Again, 


^cr — 1 _ ^cr 2^cr 

(icT + 1)^ ^ i<r + 1 ^ (Ao" + 1)“ 

=-l+e-f/* + 2(#/A)e-</Jfc (15) 

and the first two terms, when # > 3, are equivalent to 

2 (1 - «-«-«/*) - 2 [1 - -2{(t- 3)/A} e-(*-3)/*] 

= 2e-«-3)/* [1 + (2/A) (# - 3) - (16) 

This has a maximum when 


l + e-s/fc = |(^_3) (17) 

Eq[uation (17) has a root less than 5 if 

4/^ >1 + ^-2/^, (18) 

which IS an equahty if ^ = 2 7 Thus for ^ = 1 or 2 the maximum tension 
will occur before t = 5 If ^ = 1 the maximum tension is when t = S 568, 
and IS equal to 3 266 gmjw — 1 633 times the statical tension If ^ = 2 
the corresponding results are = 4 368, 2 520 gml^, and 1 680 times the 
statical tension 

The third term enters at time 5, and afterwards is equal to 
2 [1 - 2 {(t ~ 5)/kY 

If ^ = 4, the maximum stress is when ^ = 6 183, and is equal to 2 29 gmj'm 
The statical tension is 1 25 gml'nr, so that the ratio is 1 83 
Love proceeds by a method of continuation, but the present method 
IS much more direct, and probably less troublesome in application 
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4 6. A general proof that the results given by the operational method, 
when applied to the vibrations of continuous systems, are actually 
correct, has not yet been constructed It would be necessary to show 
that the solution actually satisfies the differential eq[uation and that it 
gives the correct initial values of the displacement and the velocity at 
all points The proof that it satisfies the initial conditions would be 
difficult in the most general case We have seen that in one of the 
simplest problems, that of the uniform string with the ends fixed, the 
verification that the solution is valid for the most general initial dis- 
placement IS equivalent to Fourier's theorem But for more specific 
problems the operational solution is equivalent to a single integral, and 
the direct verification that it satisfies the initial conditions is usually 
fairly easy To show that it satisfies the differential equation, it would 
be natural to differentiate under the integral sign and substitute in the 
equation But in practice it is usually found that the integrand, near 
the ends of the imaginary axis, is only small like some low power of 
1 /k (the second in the problems of 4 3, 4 4, and 4 5), and consequently 
the integrals found by differentiating twice under the integral sign do 
not converge But we can proceed as follows If a part of our solution is 


y =/(a) = A jV(^) exp .C ^ , 


and this integral is intelligible for all values of x and t within a certain 
range, we have 


g-^ = Lim 


ytJrlh + Vt- 


T I f ^ A 6'^^“ 2 + dK 


Also 


0^2 


iLim^^ 


'^yx±y^ 




= Lim 


27rt 



exp K 



W K 


The two integrals are identical before we proceed to the limit, and 
therefore their limiting values are the same Hence 




and the differential equation is satisfied This argument does not 
assume that the derivatives of y are expressible as convergent integrals, 
but only that they exist 
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The argument breaks down at points where the second derivatives do 
not exist as for instance in 4 3 at the points where the slope of the 
curve formed by the string suddenly changes At these points there is 
a discontinuity m the transverse component of the tension, so that the 
point has momentarily an infinite acceleration This is why the velocity 
changes discontinuously when a wave arrives The momentum of a 
given stretch of the string, however, varies continuously, the difficulty 
IS the fault of the representation by a differential equation, not of the 
method of solution 


CHAPTER V 

CONDTJCTION OF HEAT IN ONE DIMENSION 

5.1. The rate of transmission of heat across a surface by conduction 
IS equal to -A9 Vjdn per unit area, where V is the temperature, k a con- 
stant of the material called the thermal conductivity, and dn an element 
of the normal to the surface Hence we can show easily that m a uniform 
solid the rate of flow of heat into an element of volume dxdydz is 

k'^^V dxdydz 

But the quantity of heat required to produce a rise of temperature dV 
in unit mass is cdV, where c is the specific heat, and therefore that 
required to produce a rise dV m unit volume is pcd V, where p is the 
density Hence the equation of heat conduction is 

~{pcV) = lvW ( 1 ) 

If we put 

klpc=^h\ ( 2 ) 

h" IS called the thermometric conductivity, and the equation becomes 

dV 

^ = (3) 

In addition, there may be some other source of heat If this would 
suffice to raise the temperature by P degrees per second if it stayed 
where it was generated, a term P must be added to the right of (3) 

It IS usually convenient to write o- for djdt, and 

<t = AY (4) 

The operational solutions are then functions of g' , but o' must be expresse d 
again in terms of <r before interpreting 
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6.2. Oonsidor liint a utiifonn rod, with ite aides thermally insulated, 
and initially at teiujieiaturo H At time zeio the end $!= 0 is cooled to 
tomtieiatiiro zero, and afterwanla nmmtaiiiod at that temperature The 
end .1 - ^ la kept at teuiporatiiro iV B'lnd the vaiiatiou of tompeiaturo at 
othei pointei of the rod 

The problem being oue-dimonaioiial, the wpiatiou of heat conduction 


la 


dV 


- 0 , 


( 1 ) 


while at time 0, T is wpial hi aS' Hence the siibauliaiy equation is 

-tas; . ( 2 ) 

or • • (3) 


The end couditious are that 

r 0 when ^ - 0, 1 

V AS'whotiir -I ] " ^ ^ 

'>“■» <' “'Xi ■>)■ ■ ® 

The iiiteKnuul in an nvnn fuia^iion of and iliorofor© a single-valued 
function of <r It Im polos whcni 

ql ± Oifw ; that IH, <r *(6) 

wlioro n la any Integer But the negative values of q give the same 
valm^m of <r as the pcmitive values, and tlierofoie when we apply the 
partial«frac*iion rale we need eonsidor only the positive values The part 
arising from o* 0 is 

«('-'7)-''a 

The general term is 

"" ( - (anr) 

nrr I 

and the complete selutiou is 

r -asI'^i a ^ wl . . .. (9) 

J. n trt’f ‘ J 

If ’Tfht^ll is moderate, this senes oouvorgos rapidly, and no more con- 
venient solution could be desired It evidently tends in tbo limit to the 
steady value h'x/l. 
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But if 7rkt^/l IB small the ooiivorgcuum will \m nltnv. In thin casn we 
may adopt a tortti of the exparmion methcHl appliinl io Wo cm 

write (5) m the form 

-*0(1 "+ ^ ^ )1 . .(10) 

For if wo iiiter|)rot this as an uitoKral alcmg f lu^ patih tlu' ar^iumnit. 
of q 18 betwoou ± Jir at all point^^ (»f tlu^ path, ami the series cmwig^B 
iimformly Integration U‘rm by term is therefoie juHtifmbli\ ami we may 
interpret term by tciin Now 

qr xJ/fi ( H ) 

1 Hlf (12) 

thV 


and by 2,4 (16) 
Hence 


When w ih great, I -Erf w is wiiall ooniimn'd with c It tlieii ufj/zt^ 


is moderate, but Z/2A<^ largo, this Hori<*H is nmiiily (Mmvc'rgi'ut, »U(d can 
in most oases bo reduced to its first form 'I’liia solution is thereforo 
conToment in those cases whore (9) is not 


6,3. One-dlmandonal /l(»o of imt in » ri’ffhni infinih’ in Indh direc- 
tions First suppose that at tune 0 the distrilmtion ttf teniperaturo is 
given by 


V mr). 


.( 1 ) 


We have just seen that the function expressed in operational form by 


«■«* - I - Erf 

Vi/iti 

satisfies the differential ocpiatiou 


..( 2 ) 


3'/ 

dt 





(2) 


for positive values of at; and by symmetry it. will also satisfy it for 
negative values of a, since the function and its derivatives with regard 
to X are continuous when !»■■•() Also when t tends to zi*m this funetion 
tends to zero for all positive values of at, and ts> 2 for all negative 
values It follows from those facts that the function 


K2-<»-*”)'i[l + Brf 




•ft) 


satisfies the differential equation for all valuos of x and all positive 
values of t, and when t tends to zero the funetion tends tn zero for 


* Heaviside, MUctromagnettc Thinryt fl9 70^ 11147 
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negative values of x, and to unity for positive values Hence this function 
IS the solution recLuired 

Suppose now that the initial distnbution of temperature is 

(5) 

where /(x) is any known function Then this is equivalent to 

F=/" /(i)ci{B(i-x)} (6) 

But the solution for positive values of the time that reduces to 
when t==0 IS, 

l+Erf^1 

Hence by the principle of the superposibility of solutions the solution 
of the more general problem is 




(7) 


= 4=r /(f)e-(f-*)V^ ^ 

(8) 

Put now 

i=x + 2kt^k 

(9) 

Then 

F= 4 = r dx 

VTT 4— 00 

(10) 


This IS the general solution obtained by Fourier 


5.4. If the temperature is kept constant at x=Oy but the initial 
temperature is/(^) for positive values of x, we may proceed as follows 
If we consider instead a system infinite in both directions, but with the 
initial temperature specified for negative values of x so that 

/(-^)=-/(^), (1) 
we see that the temperature at ^ = 0 will be zero at all later instants and 
the solution of this problem will fit the actual one Hence 

F= 4^ f 7(f) Ji ±r /(I) (2) 

( 3 ) 
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If in particular the initial temperature is everjTvhere and 


V=8 



e-^'dk 


= S'Erf~ (4) 

2hti 

Thus the solution of 5 2 is regenerated In Kelvin’s solution of the 
problem of the coohng of the earth, 5 3 (10) was adapted to a semi- 
infinite region in this way 

The value of dVjdx at the end ^ = 0 is obtainable by differentiating 
the solution valid for a semi-infinite region In Kelvin’s problem 





The same result is found by differentiating (4) 

W e notice the curious fact that although the original exact solution 
for a finite region in 5 2 (5) is a single-valued function of <r, a square root 
of t appears in the approximate solution for a greatly extended region 
The reason can be seen by differentiating the exact solution It gives 


which IS a single-valued function of or But when we use Bromwich’s 
interpretation we find that 


(Sx=o = ^ Lt (t) 

which IS again single valued, but if Kk is specified to be real and positive 
when K IS real and positive, it has a positive real part at all points on 
L, and therefore if I is great coth is practically unity The integral 
IS therefore equivalent to 


fi, 

which IS our mterpretation of a-^SIh, and is equal to 8lhJ{irt) We could 
have started by specifying the sign of to be negative when k is positive, 
but then and coth kHIIi would both have been simply reversed in sign, 
and the same answer would have been obtained 


5.5. Imperfect coohng at the free end of a one-dimensional region 
"With the initial conditions of 5 2, let us suppose that the end = Z is 
maintained at temperature S as before, but that the end 0 is not 
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effectively cooled to temperature zero Instead we suppose that it radiates 
away heat at a rate proportional to its temperature At the same time 
heat is conducted to the end at a rate hdVjdx per unit area These 
effects must balance if the temperature at the surface is to vary con- 
tinuously, so that instead of having F= 0 at the end as before we shall 
have a relation of the form 

9F 

__aF=0 at (1) 

The equation 5 2 (3) is unaltered, and the operational solution is 

F= /S{1 - A sinh q {l-x)} (2) 

where A is to be determined to satisfy (1) Hence 

qA cosh - a (1 - J. sinh ql) = 0 (3) 

and r = (4) 

I gcosh^^+asmh^^J ^ ^ 

The roots in o- are real and negative, and we can proceed to an inter- 
pretation by the partial-fraction rule as usual Or, using the expansion 
in waves,” we have 


I 




If the length is great enough to make the terms involving inap- 
preciable, we can reduce this to its first two terms, thus 


L q + a J 


If a IS great, the solutions reduce to those of 5 2 , this is to be expected, 
for (1) then implies that F= 0 when ^2? = 0, which is the boundary con- 
dition adopted in 5 2 If a is small, F reduces to S , the reason is that 
this implies that there is no loss of heat from the end, and therefore the 
temperature does not change anywhere For intermediate values we may 
proceed as follows If 


Bromwich’s rule gives 


^ f ah 

'~2mJzKi + a/ 


exp ( Kt 


K^X\ (Jk 
h ) K 


Put 




( 9 ) 



60 


CONDUCTION OP HEAT IN ONE DIMENSION 


The path of integration for A, is a curve going from to Be^^% 

where B is great, passing on the way within a finite distance of the 
origin on the positive side Denote this path by M Then 




(10) 



(11) 

But 







= 1-Erf-^ 

(12) 


The second part of (11) can be written as an integral with regard to ju-, 
where 

lx = X + ah (13) 

along a path obtained by displacing M through a distance ah, but the 
integrand is regular between these paths and the route M may still be 
used The second part is therefore 



- i exp ^2aAt + exp + aa;) 



= - exp (a^A^f + aa;) |l — Erf 1 

1 2At^ j 

(14) 

and 

1 -Et£ i - exp (■/ + cue) {1 - Erf (i + y)} 

(15) 

where 

^ = xj^A'^ , y = oAt^ 

(16) 

Hence 

y=S[Etfi + exp (y=+ ax) {1 - Erf + y)}] 

(17) 

This is the same as Eiemann’s solution* 


The temperature at the end is 



Vx=o = exp y® {1 - Erf y} 

(18) 


whence the temperature gradient at the end follows by (1) For small 
values of the time the temperature at the free end falls continuously , 
the temperature gradient there is not instantaneously infinite as in 
5 4 (5) For great values of the time we can use the approximation 2 4 
(21), giving 

* Biemann Weber, Fartielle Differenttalgletchungen, 2, 95-98, 1912 
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This IK (ximvalont to ouo found by Heaviside* Heaviside gives also an 
expansion in a convergent senes, suitable for small values of t, but it 
IS probably loss convoniout than tlie equivalent expression (18) in finite 
terms W o see from ( 1 9) that the longer the time taken the better is 
the approximation to (d « given by the simple theory of 5 2 

6.6. A long rod is fastened at the onda; = 0, the other end not being in 
contact with a conductor Initially it is at temperature 0, but at time 0 
the clamped end is raised to temperature 8 and kept there Each part of 
the rod loses heat by radiation at a late proportional to its temperature 
The dilferontial equation is now 

a r as V 

(X) 


where a ib a constant Lot us put 

cr + ^^2^ 

and write the ecniation 

3^^-^r-o .(3) 

At w I there is no comlnotion out of the rod, and therefore 3F/9a? 
vanishes Also V - 8 when x 0 The solution is 

(*) 

This can bo expanded in poweis of The first term, which is the 
only one wo require, is 

F -iVir’™ . ( 5 ) 


8 [ ( (k + a'-‘)^iB') cIk 


exp 




But if 

the toim in l/(A-a) bocomos 




^ OMI (i- «-<)} oxp (- “) 


lUectTomagnetu Theory^ 2, 15 
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as m 5 5 (12) The complete solution is 


F=4*S' 



hJ I >2,hth ) 


When a ~ 0 this reduces to 



so that the disturbance of temperatme sproads along the lod, the I, line 
needed to produce a given rise of temperature at diskuieo .r fiom the 

end being proportional to ^ lit is small enough to make Huiall, tho 
error functions will be piactically unity except wheie it is not groat 

compared with For such values of .r the exponentialH are nc^arly 
unity* Thus at first the heating proceeds almost as in the absence of 
radiation from the sides of the rod 

But if is great and xj^ht^ small or moderate, the first error func- 
tion in (11) IS practically - 1, and the second + 1 Thus in these eon** 
ditions 

This IS seen by a return to the original ocpiation to be the solutiom 
corresponding to a steady state This will hold ho long as 

IS large and positive, even if xj^ht^ ih itself largo The region where f ho 
approximation (13) is valid therefore spreads along the rod with velocity 

2aA, if has once become'Jargo, each point on ilm real reaches a nearly 
steady temperature at a time rather greater than x/2ak 

* If at^ 18 large and x/2ht still greater, the solution is nearly 

F=i/Se-«»/4l-Erf®” (14) 


If further the argument of the error function, is largo, wo can write 


V = S e exp 


(- 


(fl;-2aA0“) 
U'‘t [ 


2/tt^ 
w -■ 2cLht 



-a“ij 


— 2alit 


... (I.'.) 
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In the regions that have not reached their steady state, the temperature 
resembles that for the problem without radiation, except that the small 

factor must be introduced 


5,7. TAe coohng of the earth Cooling in the earth since it first became 
solid has probably not had time to become appreciable except within a 
layer whose thickness is small compared with the radius It is therefore 
legitimate to neglect the effects of curvature and treat the problem as 
one-dimensional Radiation from the outer surface must have soon 
reduced the temperature there to that maintained by solar radiation, so 
that we may suppose the surface temperature to be constant and adopt 
it for our zero of temperature The chief difference from the problem of 
5 2 IS that we must allow for the heating effect of radioactivity in the 
outer layers Suppose first that the quantity P defined in 5 1 is equal 
to a constant A down to a depth and zero below that depth Take 
the initial temperature to be 8+ mx, where w is a constant Then the 
subsidiary equation is 

and = - ^2 ^ 

and the solutions are 


0<x<H 
H<x, I 


( 1 ) 


and 


A 

V= -^ + 8 + mx + Be-^^ + Ce^^ 

(fnr 

= 


0<x<H 

H<x 


( 2 ) 


A term in e^ is not required in the solution for great depths, because 
it would imply that the temperature there suddenly dropped a finite 
amount at all depths, however great, in consequence of a disturbance 
near the surface The conditions to determine B, O', and D are that V 
vanishes at x = Q, and that V and dYjdx are continuous at ^ = JST 
Hence 

+ 0 + /S^4*A/cr=0, 

Be-^^ + A/(r=^ 

Solving and substituting in (2), we find 


( 4 ) 

(5) 


Y=8(l — e mx + {1 — 6""^^ wAiqx} 0<x<B. 

A 

V - 8 (1- + mx + {do^qH -1) x>H 
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These solutions involve terms of the form ^ whore a is positive To 
interpret, we write 

i^a\ (!k 


1 If / ^ K^fA da 




Wi 

.1 

27rt 


\dX 


on integrating “by parts But 

^ / exp^: 
nJM \ 

and a further integration by parts gives 

* ; /h“*p ' x) V - -L L(^‘' i>) (*■’* ■ T) 






( 7 ) 

(H) 


2'n’i . 






2h 




. ( 9 ) 


InaU .CIO) 


Thus we can obtain a solution in finite terms, whuh is easily seen to 
be identical with that previously obtained by a more laborious inoth<«l *. 

The temperature gradient at the surface is found most easily from 
(6), we have 



r / K^JI\dK 


But 


1 _„rr h 

q 2vi . 


■(V2) 


A 








= 2^ (#/w)4 e-mw -u(i- Brf 


(I a) 


and therefore 
/0F\ 8 


/% 


— 4 . ^ 4- 

aCttO^ ^ 


[ir- //Erf ^ ^ + 2 A Q* (1 - 


( 14 ) 


This IS identical with the result given in earlier works of mine excopt 
that in these the factor 2 was omitted from the last tomi 

* JejBfreys, Phil Mag 82, 575-591, 1916, The Martha 84 
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In the actual problem a considerable simplification arises fronot the 

fact that H is small compared with On this account we can 

expand the solutions in powers of H, and retain only the earlier terms 
Thus for the surface temperature gradient we have 


^8 AH/, E \ 

and for the temperature at depths greater than H 


(15) 


r=Sil-e-^^) + mx + ^ 


- V- - / ■ 

= ^^ + (^_^)Erf- 


-qx 


AH^ 

2^“ 


(16) 


6 -71. An alternative possibility is that the radioactive generation of 
heat, instead of being confined entirely to a uniform surface layer, may 
decrease exponentially with depth In this case the subsidiary equation m 

( 1 ) 


at all depths We already know the part of the solution contributed 
hj S+mx The remainder is 


A 

¥ q^-a? 


( 2 ) 


But 

and 


(jf — a-^) 0- — a? W a- 




' or ^ + a Q-olJ 


C?5) 


= J [e-a* - exp (/ + a^r) {1 - Erf (X + y)}] 

+ i - exp (y"* - a^) {1 - Erf (X - y)}] . . (4) 

where 'k=xliht^ , y = aht^ .. (5) 

A 

Hence W = {exp (/ - ax) - exp (- ax)} 

+ [(l - Erf X) - 1 exp (y“ + aa;) {1 - Erf (X + y)} 

- i exp (y’* - oir) { 1 - Erf (X - y)}] 

= [1 - Erf X - 6-““= - J exp (/ + M!) {1 - Erf (X + y}} 

+ 1 exp (y’* - our) {1 - Erf (y - X)}] (6) 
which 18 the same as the solution obtained by Ingersoll and Zobel*. 


* Mathematical Tlieoi y of Heat Gonduction, Ginn, 1913 


J 


5 



66 PEOBLBMS WITH SPHEBIOAL OB OYUNDBIOAB STMMETBY 


The contribution of radioactivity to tlio tomporaturo gradient at the 
surface is 

/ 9TF\ ^ 1 . (7) 

V q-¥ a 

= /^{l-oxpyHl-Eify)}. .(H) 

When 7 18 great, as it actually is, we have 

(x—^ U ^ A- ^ V (SO 

\ 9^ /o3«o \ y ijir) k?a\ ahfJiTtty 

6.8. The justification of the method is oaaioi in piohlcmH of heiit 
conduction than in those of the last chapter, bocauBC tlie nitegraiKlH 

always contain a factor exp tends to zero when k tends 

to c± t oo in such a way that the integrals obtained by diiloreniiating 
under the integral sign always converge, and can therefore be sub- 
stituted m the differential equation diioctly But the integrals for the 
temperature are of the form 

and when we substitute m the equation 


dt ^ " 


the integrand vanishes identically 


OHAPTEE VI 


PEOBLBMS WITH SPHERICAL OR CYLINDRICAL 
SYMMETRY 


6.1. So far we have treated only problems of wave tranHimsHion or 
conduction of heat in one dimension If our system has spherical 
symmetry, the equation of transmission of sound takes the form 

= 0 , . (1) 


where ^ is the velocity potential, and 


13/ 

5 ^ — I r 

r® dr \ 


2 a4> 
0r® t Sr 


1 - 


. .( 2 ) 
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The differential equation is therefore equivalent to 

= 0 . . (3) 


This IS of the same form as the equation of transmission of sound in 
one dimension, taking the place of but differences of treatment 
arise from differences in the boundary conditions A similar transforma- 
tion can of course be applied to the equation of heat conduction 
If the symmetry is cylindrical, we take Cartesian coordinates cc, y, 
and put 

= y = xtaxL<l>i . (4) 

and if ^ IS independent of z and < 5 i» we have 




1 A 

w dw‘ 



(5) 


which IS capable of no simple transformation analogous to that just 
given for the case of spherical symmetry This fact gives rise to striking 
differences between the phenomena of wave motion in two and three 
dimensions 


6.2. Consider first a spherical region of high pressure, surrounded 
by an infinitely extended region of uniform pressure, the boundary 
between them is solid, and the whole is at rest Suddenly the boundary 
IS annihilated Find the subsequent motion The problem is that of an 
explosion wave We suppose the motion small enough to permit the 
neglect of squares of the displacements At all points 

a, 

Initially there is no motion, so that ^ is constant, and may be taken 
as zero The pressure P is -pd^/dt, where p is the density, this is a 
positive constant when r < a, but zero when r>a Then we can take 
the subsidiary equations to be 



\ 

(2) 

=0 9 >a J 

1 


The pressure must remain finite at the centre Hence 


~ + J. sinh ^ r<a 
cr c 

(3) 

= B exp (— crrjc) r> a 

(4) 


There can be no term in exp (crr/c) in (4) since it would correspond to 
a wave travelling inwards Now the pressure and the radial velocity 


5-2 
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must be continuous at r = hence ^ and d^/dr must be continuous 
These give 

-^ + Asmh^ = Sexp^-^^, ( 5 ) 


cr ^ (ra (T ( (raj\ 

' + - A cosh — = — B exp , 

' c c c ^ V c y ’ 


whence 

A = ^(e + acr) ]g =. ^ _ ^ 0 ,) gcra/c „ + ao) (7 ) 

Thus outside the original sphere 

^ (c - ao-) e-'^ (»•-«)/« - ^ (c + ao-) e-“' (8) 

The associated pressure change us 




^-<r(r~a)/c_ (r+a)/i>J (9) 


But ~ (3^ = 0 when t is negative 

= ct — a when t is positive 
Hence ^ 0 -^{r{r-a)ic == q ^bien ot<T — a 

and =ct-(r — a) — a 

= ct -7 'when ct>r-a (11) 

Similarly = o when ct<r + a 

and =ct -7 'when ct>r-\'a (12) 

Hence the pressure disturbance is zero up to time (r — a)jc, when the 
first wave from the compressed region arrives, and after (r + a)lc, when 

the last wave passes At intermediate times it is equal to ^ (r ~ ct). 

Thus it IS equal to pa/ 2 r when the first wave comes, — pa/2r when the 
last leaves, and varies hnearly with the time in between The compres- 
sion in front of the shock is associated with an equal rarefaction in the 
rear* 

Within the sphere the pressure is 
7 ^ (c + a<^) siDh ^ e-Wo] 

* Of Stokes, Phtl Mag 34, 1849, 52 
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This IS equal to p up to time (a-r)lc, then drops suddenly to 
p (1 -a/2^), decreases linearly with the time till it reaches -paj^r at 
time {a + r)/6*, and then rises suddenly to zeio The infinity m the 
pressure at the centre is only instantaneous, for the time the disturbance 
lasts at a given plaeo is 2r/c, which vanishes at the centre It is due to 
the smiultaneous amval of the waves from all points on the surface of 
the sphere, at other points the waves from different parts of the surface 
arrive at different times, giving a finite disturbance of pressure over a 
finite interval If r < the pressure becomes negative immediately on 
the airival of the distuibaiice 

The behaviour of the velocity at distant points is similar to that of 
the pressure If u is the radial velocity, 

u = d^ldr 

= (14) 

If r IS great the first term is proportional to 1/r, the second to 1/r^ 
The first is therefore the more important But the first term is simply 
a multiple of the pressure Hence there is no motion at a point until 
time (r-a)/c, when the matter suddenly begins to move out with 
velocity a/2r Tins velocity deoieases linearly until time (? + a)/c, when 
it IS -a/2r, and then suddenly ceases The total outward displacement 
contributed is zeio The second term, however, gives a small velocity 
which vanishes at the beginning and end of the shock, and reaches a 
positive maximum at time r/c It produces a residual displacement, of 
Older ajr times the greatest given by the first term, this represents 
the fact that the matter originally compressed expands till it reaches 
normal pressure, and the surrounding matter moves outwards to make 
room for it 

6,21. Consider next the analogous problem with cylindrical sym- 
metry With analogous initial conditions, the subsidiary equation is 


f 

^ dm) 

1 

11 

1 


= 0 

•nT>a J 


The solutions are Bessel functions of imaginary argument and order 
zero, and Ko(<rmlc) The latter is inadmissible within the 

cylinder, because it is infinite when = 0 The former cannot occur 
outside it, for the following reason The interpretation is to be an 
integral along a route through values of the variable with positive real 
parts, and when m is great the asymptotic expansion of lo c) con- 
tains exp (K-nr/c) as a factor Hence the solution would involve exp (o-w/c) 
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and therefore a disturbance travelling inwards The solutions are therefore 




)<'m<a\ 



^ \c J 1 

(2) 

Also d^jdt and d^jdw must be continuous at ct* = a Hence 




(3) 


-sir.' (2) 

(4) 

Also we 

haTe the identity 




(5) 

Hence we find for points outside ‘m = a 




(6) 

But* 

lo (z) = ^ exp {z cos 0) dO, 

(7) 


rco 

B'o(z)= / exp (-« cosh D)<fo, 

(8) 

whence 

Jo 


4> = - 

"271^1(5 in/o Jo ” + “ cos cosh cos 



(9) 


Performing first the integration with regard to k, we obtain a function 
of the form H (f — 6). Thus 

where the range of integration is restricted by the condition that 

+ a cos ^ - OT cosh V > 0 ( 11 ) 

It follows at once that there is no movement at a place until time 
(-zzr — a)jc Integrating next with regard to we see that the admissible 
values range from 0 to cosh“^{(c^ + cos provided the quantity 
in the parentheses is greater than unity Hence 

( 12 ) 

The notation is that of Watson’s Bessel Functions 
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So long as ct>'m — a, there will be a finite range of integration with 
regard to 0^ the lower limit is then always zero If ct>vT + a, the 
inequality is satisfied for all values of 6 up to tt, and therefore tt is the 
upper limit If the inequality is not satisfied when 0 = tt, and 

the upper limit is cos””^ disturbance at any point may 

therefore be divided into three stages, the first until cif = 'zar-a, the 
second from then till + a, and the third later 
We are concerned chiefly with the pressure This remains constant 
until Gt^m -a, m the next stage it is equal to 

jp_pa co s 0 dO 

h {{ct^a cos 0)^ - OT® 

and in the last to a similar integral with the upper limit replaced by x 
By applying the transformation 

Cif + (2 - 'zzr - 26 , + a cos ^ - 'zzr = 26 COS® if/ (14) 

and integrating on the supposition that b is small, we find that soon 
after the arrival of the wave 




When the wave arrives the pressure therefore jumps to and 

then falls by 3 cja of itself per unit time The corresponding fraction m 
the spherical problem is cja At time w/c the pressure is still positive , 
but when c# - w + « or more and 6 is greater than lir, the integrand in 
(13) 18 numerically greater than for the supplementary value of 6, and 
thus I* 18 negative The passage of the wave of rarefaction is therefore 
indelinitely piotraetod To find out how it dies down with the time let 
us suppose that ct is greater than w, and that a is small compared with 
either Then 

P = -i/> . ( 16 ) 

approximately The lesidual disturbance falls off like 


6.22. If the motion was ono-dimonsional, as for instance if the original 
excess of pressure was confined to a length 2a of a tube, the resulting 
disturbance of pressure would consist of two waves, each with an excess 
of pressure oq.ual to \p, travelling out in opposite directions with 
velocity c Comparing the results for the three cases, we see that the 
first disturbance at a given point outside the region originally dis- 
turbed, in each case at the same distance from the nearest point of it. 
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occuis in each case at the same moment The increase of proHsiire in 
the one-dimensional problem ns Ip, in the two-dimensional one 4/) («/«)* 
and in the three-dimensional one 3p(rt/ii) In the first case the iirmsiiro 
remains constant for time 2a/c, and then drops to sroto ami remains 
there In the cylindrical problem it beRiim to tall instantly, and becomes 
negative in an interval less tlwin 2(tlc, it tlieii roaclies a negative maxi- 
mum, and dies down ap^aiii asymptotically to siero In the spherical one it 
decreases linearly with the time ami roaolios a lucrative value ct[aal to 
the original positive one at tune 2rt/c, thou it suddenly hocoines zoio 
again 


6.3. Dvoerging warns produced by a iplwre osc.illating ladmUt/* 
Suppose that a sphere of radios a begins at time <) to oscillate radially 
in period 2ir/« We reiiuiio the motion of the air outside it 
The velocity potential 4> satisfies the oiinakou 

0*^ 0^ 

• ..( 1 ) 

Initially all is at rest, the solution is thoroforo 

H> Ae-"l< (2) 

When r = a the outward displacomout is, say, sm nt when t> 0, and t,ho 
outward velocity cos nt Iloiice 


cr® + 


and 


0.2 ^2 f 

(V ■+ iC^) ( 1 + crajc) ~ 0 


(;n 




(far 

c“ -s afri 


fid 

c 

na 


exp 


(-:)} 


-oxp 



r r - a) 

SlUM 

('■ .} 

J Ct 


n- 



.(4) 


when ct>>i -a 

The solution has a periodio part with a period equal to that of the 
given disturbance, together with a part dying down with the time at a 
rate independent of n, but involving the size of the sphere As there 
IS no corresponding term m the problem of G 2 wo may regard it as 

Love, Froc Land Math, Soc, (2) 2 1904, 88, Bromwich, ib (2) 15, 1916, 481. 
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a result of the constraint introduced by the presence of the rigid sphere 
Its effect on the velocity or the pressure is to that of the second term 
m a ratio comparable with (cfnay 

6.4. Aspherical thermometer bulb is imtiallyat a uniform temperature 
eq^ual to that of its surroundings The temperature of the air decreases 
with height, and the thermometer is carried upwards at such a rate that 
the temperature at the outside of the glass vanes linearly with the time 
Find how the mean temperature of the mercury varies^ 

The temperature within the bulb satisfies the equation 


^.(rV)-q^rV=^(), ( 1 ) 

where ^ (^) 

That at the outer surface of the glass is Gt, where is a constant 
But the glass has only a finite conductivity, so that the surface condition 
at the outside of the mercury is 

^ = (3) 

jK" being another constant The solution of (1) is 

F=ysmhgr (4) 

j ^ 1 _ Ct JS^Gjcr 

and ( ) gives - ^ V j 

where a is the inner radius of the glass 
The mean temperature within the bulb is 

F„ = -, V&r 


■- (S'" cosh qa - smh qa) 


_ %KG qa cosh qa - smh qa 

” aaq^ Ka smh qa + qa cosh qa - smh qa 
In applymg the partial-fraction rule, we notice that near 

Tr 3irg + 

" aaq^ qa{Ka + ^Kq^a? + ^q^a^) 




* Bromwioli, PMZ Mag 37, 1919, 407-419, A E McLeod, PM. 37,1919, 134 
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SO that there is a constant lag in the temperature of the mercury in 
comparison with that of the air 

The other zeros of the denominator give exponential contributions, 
which are evaluated in Bromwich’s paper 


6.5 ■ A cylinder of internal radius a can rotate freely about its axis 
It IS filled with viscous liquid, ai;d the whole is rotating as if solid with 
angular velocity wq The cylinder is instantaneously brought to rest at 
time t = 0, and immediately released Find the angular velocity later. 
(Math Trip Schedule B, 1926 ) 

The motion is two-dimensional, and there is a stream-function xj/ 
satisfying the equation 


vv. 




( 1 ) 


where v is the kinematic viscosity Since the motion is symmetrical 
about an axis the right side is identically zero Put 


| = <r, = = (2) 

Initially i/^ = V^ij/ = 2o>o (3) 

and the subsidiary equation is 


(V2_^) VY = -2j*’=a)o (4) 

The solution is 


’Z'- J./o(njr)+5Z'«(»-CT-) + C'logOT + Z>+ . ( 5 ) 

^e velocity must be finite on the axis, hence 5 and 0 are zero Also 
D IS independent of w and therefore cannot affect the motion 
If /be the moment of inertia of the cyhnder per unit length and a, 
its angular velocity, the equation of motion of the cylinder is 

y.d(l> 

= (6) 

where is the shearing stress in the fluid Now 

^ \daf dfj 

evaluated on the outer boundary at the point (a, 0) Hence 
A 

Ps = vp—i [ralo (ra) - 21^' (ra)] . (g) 
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Since the cylinder starts from rest, 

/o-co = - 27rvp Ar [ra Jo Q a) - 2lo {ra)"] (9) 

Also aw must be eq[ual to the velocity of the fluid Hence 

ao) = Ar Jo' (ra) + awo (19) 

Eliminating A we have 

0 ) [( JTo- 2) /o' (ra) + 7a lo (ra)] = coo [ra Jo (ra) - 2 Jo' (ra)] (1 1) 

where ^^rvpa K=I • (1^) 

The operational solution is therefore 

?aJo (ya)--’2/o'(^a) /.qn 

(K(t ~ 2) Jo' (ra) + ra U {ra) 

But 

Jo (m) = 1 + + , /o' (ra) = ira + iV + (1^) 

The contribution from or =0 is found to be 


CDo^^ 

a^ + vK 


= 0)1 


(15) 


say This is the ultimate angular velocity The other teims arise from 
the zeros of the denominator If we write ih for r, and substitute for K 
in terms of o)i, we find that these satisfy 


|(<o„ _ 0 ,,) ^ + 2 } (la) - la /o (Aa) = 0 

The solution is then of the form o)i + The coefficients A^ are 

determinable by the usual method 


CHAPTER VII 

DISPERSION 

In the propagation of sound waves in air and of waves on strings the 
velocity of travel of waves is independent of the period In many problems 
this IS not the case, waves on water afford an important example 


7.1. Consider a layer of incompressible fluid of density p and depth 
H Take the origin in the undisturbed position of the free surface, the 
axis of upwards, and those of x and y m the horizontal plane Put 


_0_ ^ ^ 
dy^^^ “^y^ 


( 1 ) 
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The velocity potential <E> satisfioR the o<iuation 

■ (a) 

At the bottom the vortical velocity vauiahos At the free Huifaco, ho 
long as the motion is only .slightly distuihod from rest, wo have 


04> Sf 
fz dt ’ 

whore ^ is the elevation of the Iroo Hurfaco ^J'hoii wo iniist have 


(3) 


(rco.Hr(a: + H) ^ 
r Hill 1 1I ^ 


(0 


The pre.ssuro just under the free Hiufaco is -7>r*C whoro 7V) is tho 
surface tension 

But by Bernoulli’s equation it is also equal to ~p(r4>- (jf>l+ 
and F(f) does not affect the motion and can ihoroforo bo omitted Wo 
have therefore the further surface condition 


-<r^ (jj~Tr')l: . . ( 5 ) 

Combining this with (4) wo have tho differential oijuation for ^ 

{<r“-(p'~7’r“)rtanrZr}^: -0. .. . (fi) 

If the fluid starts from rest with ^ equal to i'g, a known funotion of x and 
y, the subsidiary equation will have a term a-\ on tho right, and tho 
operational solution will bo 

-(y- 7k“)» ‘W 

= cosh 7Y) tan rll]\ , . (8) 

In the corresponding problem of a nmform string in 4 1 tho ooofBciant 
of t was simply pc 

Suppose first that tho original disturbance consists of an mfmito 
eleyation along tho axis of y, with no clistiirbanco of the surface 
anywheie else Then can bo replaced by y?, and (8) is ecnii valent 
to the integral 

^ " i /i®"" ~ « tan K/y }4 dn (y ) 

The integrand has an essential singularity wherever «// is an odd 
multiple of A cannot therefore cross tho real axis withm a fhiito 

distance of the ongin, but becomes two branches exton ding to + oo above 

and below the axis 
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7.2. The method of Ueepe'^t descents The elevation of the surface is 
thus expressed in terms of integrals of the type 

hl=l oxpt{f(z)}dz, (1) 

J A 

wheroy («) is an analytic function, and t is real, positive, and independent 
of z Put, following Debye, 

f{z)^Jt + iI, ( 2 ) 

thus expressing it in real and imaginary parts If the mtegral is taken 
along an arbitrary path, the integrand will be the product of a variable 
positive factoi with one whoso absolute value is unity, but which vanes 
in arguinent more and inoio rapidly the greater t is There will evidently 
be advantages in choosing the path in such a way that the large values 
of li are concentrated in the shortest possible interval on it Now if 

z- x + iy, (3) 


we shall have ay “ ® ^ ^ ^ 

It follows that B can never be an absolute maximum But it can have 
stationary points, where 


dR _dR 
dx dy 


( 5 ) 


and wo know that those points will also be stationary points of / and 
zeros oif'(z) Those points aro usually called the ‘saddle points,’ or 
sometimes ‘cols ’ Through any saddle point it will in general be possible 
to draw two (somotimos more) curves such that A is constant along 
them In sectors between those curves A will be alternately greater and 
loss than at the saddle point itself The sectors where A is greater may 
bo called the ‘hills,’ those where it is less the ‘valleys ’ If our path of 
integration is to bo chosen so as to avoid large values of A, it must 
avoid the hills, and keep as far as possible to the valleys If then the 
complex plane is maiked out by the lines of R constant through all the 
saddle points, and A and A ho within the same valley, our path must 
never go outside this valley , but if A and B ho in different 
tho passage from one valley to another must take place through a saddle 
point In tho latter case tho value of tho integral will bo much greater 
than in the former, and thorofoie interest attaches chiefly to the case 
where tho limits of tho integral lio in different valleys 

The paths actually chosen aro specifaed rather more narrow y, e 
direction of the path at any point is chosen so that |9iJ/9sl is as great 
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aa possible If xj/ is the iiiclinatiou of tlio tangout to tlio path to the axis 
of X, we have 

a/(! , a// , fi/f 

a«- (6) 

and if this is to bo a nuuioiioal inaxiiuiiui for variations m ifi 

0 -HlUl/r/ K-OSlii." 

rix ^ 1 ) 1 / 

_ SR SI 

Sn ’ C"^) 

where dn is an oloroont of l«nf?th normal to the jiath, drawn so tliat 
when d<i is in the direction of x incrt'nsinK, dn is in tli« direction of ?/ 
increasing Hence / is constant along the path. Hnoh a path is called 
a ‘hue of steepest doacent ’ There will he one in each valh'y. In general 
the limits of the integral will not themselves ho on lint's of stt'opost 
descent, but can bo joined to them by iiatlis within tJio valleys. 

In general linos of R constant through ddh*rtnit sathllo points will 
not intersect; and there will bo only one saddle point on t'ach line of 
steepest descent. For the fonnor ovnut would imply that R has tho same 
value at two saddle points, tho latter that / has, and either of those 
events will he exceptional. It follows that as wo procoed along a line of 
steepest descent R will randy reiuih a minimum and thoii prooeixlto 
increase again For if R had a nuninium SRjSx wmihl ho soro ; hut SI/Ss 
IS zero by oonstruotion, and tlioroforo tho iioint wouhl he another saddle 
point Lines of stoojiost doscoiit usually terminate only at singulantios 
of /(*) or at infinity 

Tho path of integration once chosen, tho gi cater i is tho inoro closely 
the higher values of tho intograiul will ho eoiioontratod ahout; tho siuldlo 
points. Thus wo can obtain an approximation, which will ho hotter tho 
larger t is, by considering only tho parts of tho i«ith in those rogioiiH. 
In these conditions we can tako 

/(«) /W+i(» -c„)V"fe) (8) 

where z;, is a saddle iioint. Put 

\f" (»«) I ./I ; I a " r (yj 

Then on a hue of steepest doscent 

/(») no) 

and we can put 

arg(a:-«,) 
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on tlie side after passing througli Zq Then 


5^=2 r exp {j?/(5;J}exp(-|-J[i^r®)rf(rexp ta) 

Jo 




( 12 ) 


To get from JL to 5 it may be necessary to pass through two or more 
saddle points, with probably traverses m the valleys between the lines 
of steepest descent Then each saddle point will make its contribution 
■fco the integral 

The error involved in this approximation arises from the terms of the 
third and higher orders omitted from (8) Its accuracy therefore depends 
on exp having become small before exp (^^o)!) has 

begun to differ appreciably from unity Hence 



= 1 


i I/' (^o) I 


( 13 ) 


must be large In most cases the approximation is asymptotic, and does 
not represent the first term of a convergent senes 


7.3. In problems of wave motion we often have to evaluate integrals 
of the form 

<j6 (k) exp (k^ - rO dK, (1) 

Jjj 

•where <f> (k) and y are known functions of k As a rule (k) is an even 
function, and y an odd one When k is purely imaginary y is also purely 
imaginary We require the motion for large values of and possibly 
also of 07 

The function <f> (k) usually introduces no difficulty It does not involve 
07 or and therefore when these are large enough it can be treated as 
constant throughout the range where the integrand is appreciable 
It IS usually convenient to replace k by t/c and y by ty, and to consider 
the equivalent integral of the form 

^ — i I i// (k) exp I (ko^ - yt) d/c .. . (2) 

2tT 7 — oo 

The saddle points are given by 

a7-'yt = 0 ( 3 ) 

so that a given ratio (jcjt specifies a set of predominant values of ko But 
y' IS an even function of k, and therefore if k© is a saddle point, — kq will 
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be another, and if the adopted path passes through either it will pass 
through the other We may take kq real and positive Also 

f (^ o ) =*“^70 t ( 4 ) 

Thus if IS positive /c — /cq has argument on the line of steepest 
descent, and the contribution from /c^ is 

1 / 27 r 

^ \\W^\) ^ exp t (k„ a; - yo t) ( 5 ) 

The contnbution from - k, (where y" is negative) is similarly 

✓ 2 - 7 r 

^ ‘ ^ - rot) ( 6 ) 

and the two together give 


( 2 

^ I jq' t \) ^ (^0 ^ “ 7o ^ 

Similarly if y^' is negative the two saddle points give 
/ 2 

(^<0^-70^ + i^r) 


( 7 ) 

( 8 ) 


7 . 31 . These formulae, due to Kelvin, are the fundamental ones of the 
theory of dispersion Consider first the cosine factor, and suppose w in- 
creased by Sx and t by St Then KQX—yot is increased by 
/co8a7~7oSi^ + (i2?— yo'if) 8/Co, 

the term in Sk^ appearing because kq is defined as a function of x and t 
by ( 3 ) But the coefficient of Sk^is zero by ( 3 ) If then t is kept constant, 
t will vary with x with period 27r/fCo, and if x is kept constant, C will 
vary with t with period 27r/yo Hence 27r/fCo is the wave-length, and 
2^/70 the period, of the waves passmg a given place A phase occurring 
at a given place and time is reproduced after an interval St at a place 
such that 8 x = yo St/K^ Hence yjKf, is the velocity of travel of individual 
waves It may be denoted by c, and called the wave- velocity 
But Ko has been defined by the equation 

X yQ t = 0 

80 that a given wave-length and period always occur when m/t has a 
particular value, they seem to travel out with velocity yo', which is 
called the group-velocity It may also he denoted by C In general the 
wave-velocity and the group-velocity are unequal, so that a given wave 
changes in period and length as it progresses They are evidently con- 
nected by the relation 

S-(‘^o) = C 


( 9 ) 
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7.4. Eeturmng now to 7 1 (9) we can separate the hyperbolic function 
into two exponentials, which will represent wave-systems travelling out 
in opposite directions One of them is equivalent to 

1 

expL(KX-yt)dKj 

4:7r J-00 

(1) 

where 7^ = ^ tanh k H 

When K is small. 

(2) 


(3) 


(4) 

When K IS great, 

(5) 


In all ordinary cases Tjg is insignificant m comparison with Hence 
for some intermediate value of k the group- velocity is a minimum ^ it 
tends to infinity for very short waves, and to a finite limit for very long 
ones Three cases therefore arise If xjt is less than the minimum group- 
velocity, there will be no saddle-point on the real axis, and the disturb- 
ance will be small* If it lies between this minimum and two 

(positive) values of k will give saddle-points, and each will contribute 

to the motion If it is greater than the only saddle-point will 

correspond to a short wave The disturbance at a given point will theie- 

fore be in three stages In the first, leading up to time xl{gH)^y only 
very short capillary waves will occur Then long gravity waves will 
arrive, the wave-lengths of those reaching the point diminishing as time 
goes on Superposed on them are further capillary waves, their length 
increasing with the time At a certain moment the wave-lengths of the 
two sets become eqLual This corresponds to the arrival of the waves 
with the minimum group-velocity From then on the water is smooth 

7.41. Two typical cases therefore arise according as the wave-length is 
large or small compared with the one that gives the minimum group- 
velocity Take first gravity waves, such that xjt is small compared with 

* It can be shown that the saddle-pomts are so placed that the relevant one 
contributes an exponential with a negative index to the solution This is almost 
obvious from considerations of energy 


J 


6 



82 


DISPEESIOIT 


{gE)^ but large compared with the minimum group-velocity In these 
conditions we can write simply 

c=={glKi^, dCldK = -l{gl^)i ( 1 ) 

ihe solution is then 

( 2 !^) ^ ("g) + (2) 

f'l©* w 

and the amplitude increases towards the rear of the wave train like 

7.42. Take next the capillary waves, short enough for gravity to be 
neglected Then 

c = {TKf, G=%{TKf, dCldK = ^{TlKf ( 1 ) 

At a given instant the amplitude is therefore proportional to or to 

0 ^ The front of the disturbance is therefore composed of a series of 
oapiUary waves whose amphtude tends to infinity, and the time taken 
for them to arrive is mfimtesimal 

This impossible result arises from the form assumed for the original 
displacement In taking ^0 =i? we assumed that unit volume of liquid 
was originally released on unit length of the actual line ^ = 0 m the 
surface The mean height of this mass of liquid was therefore infinite, 
and its potential energy also infinite The system being frictionless, this 
- energy must be present somewhere in the waves existing at any instant, 
and infinite amplitudes are therefore a natural consequence of the initial 
conditions If instead we suppose that the same volume of fluid was 
originally raised, but that it was distributed uniformly between £o = + ly 
its elevation was 1/21 m this range Expressed in operational form this 
gives 

(a) + 1) - H (w - - ^"^0 ( 2 ) 

The appropriate solution can be found from 7 4 (1) by introducing a 

factor or - into the integrand If the solution already 

found makes kqI small, this additional factor will be practically unity, 
and the same solution will hold Waves whose length is large compared 
with the extent of the original disturbance will therefore not be much 
affected by its finiteness 
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But, if IS laigu wn luiwt (wuHiilor miiiarately tho contributions from 
tbo t-onus in ami 'rii.i fiic'tnr HkI will giro an extra mall 
the holufiuiiH Tim lOHult will bn that at a given instant the amplitude 

of a gravity wave will vary like ^ ami that of a capillary wave like 

WavcK whitHe lengtli ihMhert pomparodwith that of the original 
(listnihaiu'e me tlun'etuie Imavily reduced m .uuplitude. 

On dcepwatei Iheniininmmgieup-velecityis 18cm./Hoc .conesponding 
k) a wave length ef 1 1. cm and a wave- velocity of 28 cm /sec If the 
oiiginal distuiham e has a hon/,onlal extent of I cm or so, only wayos 
with lengtliH under I mii. or so will be afleetod, and the amplitudes of 
botli gravity and cnpdlary waves will incieaHO steadily with diminishing 
group-velocity A wave of huge amphtudo will theieforo bung up the 
rear, ami will leave .smooth water behind it This is obsorvahlo m the 
waves eatmed by numhops and ot hor very concentrated disturbances 
But d the extent ol the urigiiwl disturhaneo exceeds a few centimetres 
the eaptllary waves piodtieed will 1 h‘ very small, and the largest 
amplitude will he associated with a wave whoso length is comparable 
with the wnltli of tiie distiuhisl region. The largest wave pioduood by 
the splash of a brick, lor instanee, has a length of the order of a 
foot 

7.6. 'I’wo exeeptional i-ases may arise in tho troatmont of disiiorsion, 
which are both tllm.tr, ifcd in the present problem Tbo validity of tho 
approximation 7.2 (12) depends on expj/,/(i>)) being propoitioual to 
exp( * .Ur') on a line of .steepest descent has varied by a 

considerable fraction of .1 before this exponential has become small tbo 
approximation will not he good Tliis may hajipon if A is itself small, 
or if there are two .saddle- points close togetlioi. Instances occur when 
there is a waxitintm or mininiuni groiiii-veliHiity, or if tlio group-volocity 
toiuls to a finite limit when k heetunes very small In tho former case a 
value of Jilt a little greater than the minimum groiip-volocaty will give 
two slightly difierent finite values of «„ In the latter sj may bo small 
enougli for the {iroxiinity of k„ to alTeet tho contributions of both. 

7.61. Hinee y is an odd function of wo may suppose that when k is 
small 

y c»K -Cjjs’ t O^k") (l) 

and then 7.3 (2) is eipiivalent to 

I exp ‘(s,C-C„S« + CaS-'/)o?« ■ -(2) 

(>»% 
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OoiiHP(iucntly {hodiHturbiinca wa aro oousuloting produces au immediate 
nso of t.lie level of tin* water at, all dthtances greater than c^t, though 
tins rise iH wry Hinall at great distances 'I'lio maximum rise of level is 
whore ./ is rather li*ss tlian e„t, m that it has not travelled out from the 
origin so far ns a wave with the hmitiiig velocity would The maximum 
IS followed by a senes ot waves of giadually dimiinshiiig length and 
aiiipht iide, merging tilt imatelv into gravity waves of the deep-water type 

7.62. Ill t lie case where there is a minimum gioup-volocity for afinite 
wave-Iengtli, let us, with a somewhat (lillercnt notation from that used 
so far, denote the lumtmum group velocity by y,,' and the corresponding 
values of k and y by and y,. Put 

« .. (1) 

V Yo » yt's! 1 0 ! Jy.t'"si'‘ 1 . (2) 

Then 

with an analagnuH contuhiition from iho nogativo valuoB of k, honce 
with t!m Bama typn nf apiin^xmmtmu an Imloic wo nhall have 

2 r ^ 

f «/'(N,)eos(K„.r- y./}(,ly„"7} coh (v'‘ - mv) do . (4) 

•a Jo 

wliere m (5) 

yrt)^ 

The Kolutum m thorciforn the prntlnrt of a ooBino and an Airy integral, 
tlm interval hetween c*ciim(*enfivn rnmn being much greater for the 
former than for tho laitcr. In the neighhourliood of a point that has 
tmviillotl out with the minimum group-velocity the waves have the 
eorresponding pi^riod and wavelength, but thmr amplitude falls off 
rapidly towartln the u*iu\ In trout of tliis point the amplitude increases 
for a 'lUiilc and then o«ciltah»g 

In lioth tlioBeexc4»pt4onnl cuhoh we notice that the amplitude associated 
witli the critical velocity falls off only like the luvorHO cube root of the 
time, wliereas in the typiiud caMi it falls off like the inverse Bq[uare 
root* Ilenci^i the furtlicr the dmiurbance progresHOB the more will the 
waves with the oritical gtoup velociticH predominate in relation to the 
oilierH In waves on ^ater, however, these phenomena are often modified 
by the greater daiiiping effect <if viBctmity on short waves 
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CHAPTER VIII 

BESSEL FUNCTIONS 

8.1. The Bessel functions of imaginary argument are defined by the 
expansions 


2 


( 1 ) 


r=o ^ ' (n + ry 

where (n + ry is to be interpreted as r (» + »• + 1) if is not an integer 
rut 




( 2 ) 


Then 


In(2a;i) = a!-i'^ 2 




: 0 >'’ (n + ry 




~{n+r) 


r=0 »■' 

= expp-\ (3) 

where j) denotes djdoo Hence 

i3-"expp-i = a;i»7„(2a;i) (4) 

Differentiating and snbstitutmg for so we obtain the familiar recurrence 
relation 

^ In («)} = /»_i (z) (5) 

The expression on the left of (4) is equivalent to the integral 

‘ ( 6 ) 

if m IS positive, and if the path L is replaced by a loop M, passing 

around the ongin and extending to , the result holds without restric- 
tion on n Hence 

/.» = (»., (7) 

which is equivalent to Schlafli’s form* Putting also 

( 8 ) 

w.We p, 

* Watson, JBessel Functions, 175-6 
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pr(ivi(lt‘<l t.ho real part ttf ; ih piwitive If 

2 (^ ' a) ^ 

^ t (/‘‘ - 0*, . (U) 

the Ilf thi> reut lieiiij; (h'tenmned hy tho fact that A aiid /a must 
appKiach i / we find 


giviiiKi n\ ojH^mfioim! fonu^ 

/»(»r) 

IP ^ (p 

Thin eaii cwily li^s venifuiil for // 

i>f p. 


<h‘ 

0'*'- o’ {/* + (/*'““ OH"’ 

(12) 

p 

' o'-l" (/>■■' -0^ 

(13) 

0 hy expanding in negative powers 


8»2. Till* iiiloKnuul ui HI (l!2) \im bianch ponitH at ±1, but no 
oilic'r HiugularitH*H in tin* hmin partoi tin* piano The path can thoroforo 
he nnHtifn‘d to \im\H from each pahsing around ono of ihoBO 
blanch (’‘onsidor finst /„(:) The loop around i I givcH a con- 

trihutuui donotcd, tn Ucavwidea noiatnui**^, hy \ Wo have 


ILU) 



( 1 ) 


and if is grent tlie t‘xpoiieniint m apprmalilo only when It{p) ih near 
unity Then wn ran put 

M 1 e, (2) 

ninl wiiti* 




' ■' I" •' , 

^ f I'’)* 


I 

Wi 



I I 


V I 




fir, 


.(3) 


where, if the hmp Ik* Niippimed itulellidtoly tiarrcw, wo imint take (- 2e)^ 
p(mit iv(' hna'j;iiiary sin t he upper Hide sif the path, and negative iinagmary 
tin the lower, llettes* 






* TIuH^rn^ 2 , 4 l! 4 n 
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Similarly we may consider a loop around - 1 and define 

^0 («)=-[ exp {ixz) — , (5) 

the positive sign bemg taken for the root when - 1 < /i < + 1 Then we 
put 

^ - 1 - V, (0) 

and obtain 



By substituting the loop integrals for H, and in the differential 
equation for 1^, namely, 


^ 0 , 1 ^ 

dz^ z dz 


-4=0, 


( 8 ) 


we find easily that they both satisfy it and therefore constitute a pair 
of independent solutions The expansions (4) and (7) are asymptotic, 
but a fuller discussion is needed before a limit can be assigned to the 
error involved m stopping the expansions at a given term Their physical 
interest is that if the variable a is crnr/c, where w has its usual meaning 
as a cylindrical coordinate and o- denotes djdt, E, contains as a 
factor, and its mterpretation by Bromwich’s rule wiU have exp /c (# - w/c) 
as a factor in the integrand, where the real part of k is positive It 
therefore represents a diverging wave Similarly Ho represents a pure 
converging wave They are also intimately connected with the two 
Hankel functions Hi^'> and If arg » increases till z=i.y, where y is 
real and positive, the loops must be swung round till they pass to + 1 oo 
instead of- «> Then 




-r,W=i;(.s)=Q) 'j 

= -t4;P)(y) (10) 

To express 4 m terms of H and Eo, we suppose that the loop 
around + 1 passes to - co on the upper side of - 1 The integral around 
this loop IS ^Ho by definition Taking the othei loop next, we see that 
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( ft ?- iiCKativt! imaKJiHuy 1hH\v(«>u the braiicli points, and the 
int.osral is 

1 

‘iirt 

Thus /„(;) I (12) 

(JcmiparhiK with ( I) ami (7) wi> wanu to have the anomaly of a purely 
u'al function hcin^jt cqu.d to tlm muu of a pmely real one and a purely 
iniaKimiry ono Tlio «‘\jilaimtioii h that while /„ and K„ are purely real 
by definition when i w real and positive, we have had to define JIq by 



H. ClonttHun ht and 


iiieauH of a loop piiiisiiiK 1 on t he positive imaginary Hide When ju + 1 
IK iiufiginitry (/<’ 1)* has a real pait, and therefore //„ w not purely 

real. If we had defined //., by nienim of a looji paHHing - 1 on the under 
Hide we Hhould have hail to reverMo the sign of i m the oipiation 

l/M) a/„(0 ‘h\iz) (13) 

The risiKon thi.i phenmiienon doex not hIiow in the aHymplotio expansions 
w that when z w great (uiougli A'„ ia Huiallei than any term in tho 
RKyuiptotie expansion of //„, and therefore cannot atfoct this expansion. 

The ordinary Bessel function /„(c) can bo defined by 

JJz) /Jn) (14) 

i!//„(<Ot<A'„(.s)} (15) 

W“’(s)l (16) 


! (o) 


) ' coH (*” iir) t higher terms .. 
In operational form it ns given by 


( 17 ) 


(pM 1)2 


( 18 ) 
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8.3. In tho case of funotiotm of other orders, the extra factor m 
8 1 (12) TR {fj. + (/j? — 1)^} which iH unity for /t 1 and ft- Wwt I 

Hence we can define an //„ and a A';, the latter in leal fonn, and the 
former with a real aHymptotic. cxpaiiKiou. If thiH ih done the fiint toniig 
of tho expanaionH will bo uidoiiendent of w, a known lesult But when 
we proceed to the ordinaiy BohhoI function the complex faetom afloct 
the argument of tho coHino, wo have indeed StokeH’n appiovimation 

/ 2 \ i 

rOH I ir) » ‘ fl) 


8.4. HeaviBule’H procedure ih to dofiuo /4 (x) liy the opoiator 


2 ;; 


and A'„ by 



m tho present notation* Tho uiiHatiHfmitory foaturo of tboHO definitions 
is that wbothor the fuudamontrd interpietafion of an operator is an 
oxiiansion m powers of 1/p oi a coiiipl(‘\ nitegrnl, the formei operatoi 
means 2/„(ir) and the latter ±2i/,|(u‘) But Heaviside applies a process 
equivalent to changing tho vanablo from /* f,o e as alsive, and then 
expanding in asconding poweis of this now variable. 'I’bo lesiilting 
operators are equivalent to 8 2 (.1) and (7) above, the method really 
amounts to evaluating the portions of tlu* integrals tliali arise from till? 
uoighbourhoods of tho brauoli points, and tlu'refore gives //„ and A', 
correctly But the operatoi s uitroiluood at tho start do iioti represent tho 
functions Hoavisido thushca with. 


8.6. The prosont expansions can bo applied to tbo solution of tho 
problem of 6 2t. Tho opomtioual solution was 



^0 (»). as hero defined, in aooordsnoe with Hcavisltlti’s praetioo, differs from the 
deflmtion in 0 21, which is that adopted by WatBOii, following MaodonaUi, Watson’M 
/fo IB ir- times Heaviside’s. Watson oommentH on the fact that tho oatra factor 
obsoures the relation between »■„ and the Hankcl functioim, hut homtatos to remove 
It booause the ftinotion is already tabulated But iM>rhapH tables ant of hw import 
anoe than analytical eonvenutuoe, wlien Bessed functions oeour Inprobloms they 
more often than not seem to disappear from tho final answer 
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and the pressure is 


P = |p 


Q-cr{v^-a)le I 


whioli is zero till time {'m - a)lc, then jumps to |p and proceeds 

to fall by %c/a of itself per unit time 

8-6. In much of Heaviside’s work use is made of what he calls the 
generalized exponential function, defined by the series 

oo ^n+r 

S (1) 

r^-oo{n + ry ^ ^ 

where ^ is a proper fraction The series formally satisfies the differential 
equation 


and reduces to the oidmary exponential series when % is an integer 
Its mathematical peculiarity is that the part of the series corresponding 
to negative values of n is always divergent As usual {n-^ry is to be 
interpreted as r(w + r + 1) This series was recently found by Mr A E 
Ingham and myself to possess an asymptotic property Suppose we start 
with a given term, say ? = ~ m, where ^ is a positive integer, and con- 
sider the series 

8= 2 ^ (3) 

_ 1 r 's ^ 

-J_ f 

2ir6 Jt /c"-™ (« - 1) ^ ^ 

if L IS so chosen that | #< | > 1 at all points of it Then the integrand has 
a pole at k = 1, which makes a contribution ^ to the integral The 
remainder is equal to 

^ ^ ~2mj k”-’^(k-1) 

taken around a loop surrounding the origin and passing to ~ go at both 
ends This is equivalent to 


1 yjn.’^n 
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and the integral is less than 


f 


6-^^ fjL^^diJL=: £27“^”^“^+^) T(m-n + l) 


X 


n-m-i , 


r (^ — ^) sm (n — m)7r 


Hence 


\S^^\ 


so that the error caused by starting at a definite term is less than 
the preceding term, and decreases indefinitely when x becomes great 

8.7. As has been indicated, Heaviside’s work is largely concerned 
with the use of asymptotic series, which he justified mainly by appeal 
to actual computation Less defence of such series is now necessary, 
pure mathematicians have proved their validity in many cases, and 
even used them themselves Some of HeavivSide’s senes, however, are in 
their actual application convergent Thus, if we consider the asymptotic 
expansion of (x), namely 




+ (-l) 


1 3^ 

8a; 2 1 

1" 3" (2n- 


■If 


)■ 


(^Sxy^ 

this is formally divergent as it stands, but in its application to prob- 
lems of cylindrical waves the argument is o-zzr/c, and the general term 
(allowing for the factor is 


= (- 1)" 


1 1 


= (-i3« 




=(- 1 )” 


3 2n—l/ct 


Jir (2n + 1) 2”" V 2 mJ 
The general term in the interpretation of (ctstIc) is therefore 
,W^21 3 2n-l 


(- IT 


'c t — 

2w ) 


TT ?^^(2w+l)2” 

and the series converges like a binomial expansion provided ct<Z'^ 
Though a divergent senes is used in the course of the work, it leads to 
a convergent answer 

8.8. Eut the operational method raises numerous points concerning 
the relation between Pure Mathematics and Mathematical Physics in 
general Chapter I of this work sufhoes to show that the operational 



BESSEL EUNCTIONS 


93 


raetliod gives the correct results when the system considered has a finite 
number of degrees of freedom ^ but as soon as we begin to consider con- 
tinuous systems we find that the solution involves operators not 
definable in terms of definite integrations, and the method of Chapter I 
IS no longer available for a justification Bromwich’s introduction of 
complex integration serves three important purposes it enables us to 
prove in large classes of cases that the answers obtained to problems of 
continuous systems are correct, it provides a formal rule for interpreting 
operators in general, and it gives the answer, in cases where it involves 
new operators, in a form convenient for direct evaluation by contour 
integration It seems curious that physics, dealing entirely with real 
variables, should find it necessary to use the complex variable to solve 
its problems in the most convenient way The use of conformal repre- 
sentation in two-dimensional electrostatics and hydrodynamics is of 
course exceptional, and arises from the fact that the general real solution 
of Laplace’s eq^uation in two dimensions is the real part of a function of 
^ -h w/, for other equations there is no such explanation The real reason 
IS that the solutions of the linear dijGferential equations of physics are 
expressible as linear combinations of analytic functions In the method 
introduced by Bromwich himself, for instance, k is written for d/dt before 
solution, and then for the general value of k the solution is found that 
satisfies the differential equation and the terminal conditions Such a 
combination of these solutions as will satisfy the conditions when if = 0 
IS then constructed lu the form of an integral with regard to k, this 
integral is the same as would have been obtained by solving operation- 
ally and interpreting according to Bromwich’s rule The utility of the 
complex variable then rests on the fact that analytic functions of it 
satisfy Cauchy’s theorem I do not think, however, that the complex 
integral should be regarded as the definition of the operator In 2 1 (2) 
the IS the fundamental notion, and the integral merely expresses 
a convenient rule for evaluating it The rules of Chapter I m my opinion 
come nearer to the ultimate meaning of the operators, but their exten- 
sion to the operators that arise in the discussion of continuous systems 
awaits further investigation 



NOTE 


ON THE NOTATION FOE THE ERROE FUNCTION 
OR PEOBABILITY INTEGRAL 

The notation given on p 26 is one that I adopted in 1916 under the 
impression that it was in general use, and I have since used it in several pub- 
lications* My definition was recently queried by a correspondent, and I have 
not succeeded in tracing its origm I have, on the other hand, discovered a 
suipnsing confusion of other notations The earliest, due to Gausst, is 


VTT J 0 


e-^'‘ dz 


No name is given to the function hy Gauss, and there is no sign that he meant 
the notation to he permanent Of modem writers, Oarslaw, Brunt, and Ooohdge 
use this notation Fourier J gives 






but has, so far as I have traced, no modern foUowers Jahnke and Emde in 
their tables, use ’ 


f' 

Vtt J ( 




The same notation is widely used by other 
German writers Whittaker and Robinson § also use it and call the function 
the Error Function The notation 

Erfa:= [ e-^ dx 
J X 

was introduced by J W L Glaisher||, who also usedT 

Erfo x=l^ dx=lJ^ - Erf ^ 

The latter function is also called erf a; by R Pendlebury** 

* Fhil Mag. 32, 1916, 579-585 , 35, 1918, 273 , 38, 1919 718 m T<r t> a nr, 

t Werhe, 4, 9 First published 1821 
t Thione Analytique de la Chaleur, 1822, 458 
§ Calculus of Observations, 1924, 179 
il Fhil Mag (4), 42, 1871, 294-302 
IF Loc cit 421-436 
** Loc cit 437-440 



Whittaker and Watson^ use Glaisher’s notation with the meanings of Erf 

2 

and Brfc interchanged, while Jeans’s t erf^ is —= times Glaisher’s Erf 

V TT 

Numerous otlier writeis use the integrals, but omit to give any special symbol 
or use only the non-committal “ I ” 

It can hardly be denied, in view of the wide application of the error function 
to thermal conduction, the theoiy of errors, statistics, and the dynamical theory 
of gases, that it merits a distinctive notation It is equally clear that none of 
those yet used has obtained general acceptance Of them, it seems that those 
involving only the single letters 0, ^ are wholly undesiiable In dynamical 

problems these letters are in cotitinual use for couples, while ^ is often wanted 
for a velocity potential and ^ ov -xj/ for a stream function , is of course also 
often needed for an angle Further, 0 has an established meaning already in 
the theory of elliptic functions, namely, the Theta Function of J acobi, while 
^ has another meaning in the theory of numbers To avoid confusion in some 
of the applications of the function it seems necessary to use a combination of 
letters, and the only one in frequent use is Erf 
As to what function should be denoted by this symbol, the most convenient 
IS indicated by the practice of the compilers of tables, such as J ahnke and 

2 

Emde and Dale, who tabulate I du This is an odd function and 

TT J 0 

becomes unity when ^=^ 00 , two properties that make for analytical convenience 
and are connected with the fact that this form, including the factor 27r ” 
usually occurs as such in the solutions of the relevant problems Accordingly 
I think that, of the various notations proposed, 


6““*^ du 


IS the most convenient, and is worthy of wider adoption, though I seem myself 
to have adopted it under a misapprehension oiigmally 


* Modern Analysis, 1915, 335 
t Dynamical Theory of Gases, 1921, 34 



INTEEPEErATIONS OF THE PEINCIPAL OPBEATOKS 

TT/*\ rv . « 


a- ”'S{t) = 0 t<0 


=iT ^>0 




in- 




j 7 r-= “ 

(<r-o)« (»-l)t® ’ 

na 0-2 

e-<^^f{t)=f{f-x) 

fj-h f 1 s _i 

V'T?’ " =-^7' ’ " ^=2 


If <T = hY, 

e-2'*=l-Erf 


a; 1 

qe~^= — i_ 


%ht^’ " A(;rOi‘ 

1 a-2*=2A a-Vm- ^ (^i _ Erf , 

ae~~^^ X r f 

P 

{f5 + (fj"-l)4}’‘(_p2-l)i + (^2+ = W 
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